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Basic aspects of the background of gravitational waves and its mathematical characterization 
are reviewed. The spectral energy density parameter commonly used as a quantifier of the 

background, is derived for an ensemble of many identical sources emitting at different times and 
locations. For such an ensemble, 0,(f) is generalized to account for the duration of the signals and of 
the observation, so that one can distinguish the resolvable and unresolvable parts of the background. 
The unresolvable part, often called confusion noise or stochastic background, is made by signals that 
cannot be either individually identified or subtracted out of the data. To account for the resolvability 
of the background, the overlap function is introduced. This function is a generalization of the duty 
cycle, which has been commonly used in the literature, in some cases leading to incorrect results. The 
spectra produced by binary systems (stellar binaries and massive black hole binaries) are presented 
over the frequencies of all existing and planned detectors. A semi-analytical formula for il(f) is 
derived in the case of stellar binaries (containing white dwarfs, neutron stars or stellar-mass black 
holes). Besides a realistic expectation of the level of background, upper and lower limits are given, to 
account for the uncertainties in some astrophysical parameters such as binary coalescence rates. One 
interesting result concerns all current and planned ground-based detectors (including the Einstein 
Telescope). In their frequency range, the background of binaries is resolvable and only sporadically 
present. In other words, there is no stochastic background of binaries for ground-based detectors. 



I. INTRODUCTION 

The gravitational wave background [312] is formed by a 
large number of independent gravitational wave sources. 
This paper focuses on the background produced by coa- 
lescing binary systems. These are isolated pairs of mas- 
sive objects that inspiral towards each other by emitting 
gravitational radiation until they coalesce. 

We review the characterization of the background, for 
which the spectral energy density parameter, or simply 
spectral function, H(/), is often used. This function gives 
the present fractional energy density (per logarithmic fre- 
quency interval) of gravitational radiation at an observed 
frequency /. A formula for f2(/) is obtained in a clear, 
self-consistent way, for an ensemble of many identical 
sources emitting at different times and locations. This 
formula is generalized to distinguish whether the signals 
are resolvable or not, or whether they are observed con- 
tinuously or sporadically. 

The resolvability of the signals is an important topic 
of this work. Roughly speaking, signals are unresolv- 
able if their waveforms are observed simultaneously at 
similar frequencies (differing less than the frequency res- 
olution). Many unresolvable signals form an unresolvable 
background. If such a background dominates in a certain 
frequency interval, one cannot see the waveforms of its 
components, but a pattern that rather looks like instru- 
mental noise. For that reason it is often called confusion 
noise. The other part of the background is resolvable. 
The waveforms of the resolvable part can be distinguished 
and in some circumstances subtracted out of the data of 
a detector j3] H] . 
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There are many studies in the literature about astro- 
physical sources that contribute to the background at 
present. A few examples of these sources are: core- 
collapse supernovae [5] , rotating neutron stars [fH [7J , for- 
mation of neutron stars [5J [DJ , inspiralling or coalescing 
stellar binaries [1PJHT4] , inspiralling or coalescing massive 
black hole binaries [EHUD] and magnetars [T7J. But there 
are inconsistencies in the literature; for example, accord- 
ing to [15], the results of some of the previous papers [HI IS] 
(and also [T31ITD], as pointed out in [3D]) are incorrect, 
due to a wrong (1 + z) factor in the calculations. Besides, 
according to [21] , the definition of the spectral function 
used in many papers is misleading or misinterpreted. Fi- 
nally, as we discuss further on, in some of the mentioned 
papers, the continuous and unresolvable backgrounds are 
not properly defined. To avoid possible misunderstand- 
ings or mistakes we tend to present all calculations and 
definitions as clearly as possible. 

We calculate the contributions to the background of 
the strongest emitting binary systems. These are the 
ones composed of white dwarfs, neutron stars, stellar- 
mass black holes or massive black holes. The resulting 
energy spectra are given as maximum, most likely, and 
minimum expectations, taking into account the present 
uncertainties in the quantities involved. 

We show that ground-based detectors do not encounter 
any unresolvable background in their frequency window 
(the frequency range in which they reach their optimal 
sensitivity). This applies to present detectors, such as 
TAMA300 [33], GEO600 [23], Virgo [24], and LIGO [25], 
but also to planned detectors, such as the advanced ver- 
sions of LIGO and Virgo [351 US], LCGT [37], and ET 
[28] , At these frequencies, there is not even a resolv- 
able continuous background, i.e., signals are not always 
present. Whether or not these signals can be subtracted 
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out of the data is an issue we do not deal with in this 
paper. 

This paper considers a frequency range wider than the 
frequency windows of ground-based detectors. It includes 
the windows of all existing and future detectors, such as 
LISA [21] or BBO [3D], and also reaches the frequency 
range of interest for the PTA [31] . 

The obtained unresolvable background turns out to be 
dominated by white dwarf binaries (at frequencies below 
~ 0.1 Hz) and by massive black hole binaries (below ~ 
10~ 4 Hz). This confusion noise could enter the band of 
LISA and would certainly enter the band of BBO and 
the complete Parkes PTA [52"] . 

The outline of the paper is as follows: 

In Section [IT] we first explain our notation and give 
some relevant terminology. We then give a simple heuris- 
tic proof of the fact that there is neither an unresolvable 
nor a continuous background in the frequency window 
of ground-based detectors. The formula of the spectral 
function is derived for an ensemble of many identical 
sources emitting at different times and locations. The 
obtained formula agrees with that of [33]. The concept 
of resolvability is studied, and the spectral function is 
generalized to account for it. To get to that point, we 
introduce the overlap function Af(f,Af,z). This func- 
tion gives the average number of signals observed with 
redshifts smaller than z and frequencies between / and 
/ + A/, where A/ is the frequency resolution. We then 
use the overlap function to distinguish the continuous 
and discontinuous parts of the background. 



In Section |III| we present the models used to quantify 
the background of stellar binaries and massive black hole 
binaries. The main physical quantities involved in the 
calculations (such as mass ranges and coalescence rates) 
are presented in this section. A semi-analytical formula 
for the spectral function is derived in the case of stellar 
binaries. 

Section [IV] contains the main results of the paper. The 
spectral function is shown in the different regimes of re- 
solvability and continuity. The curves in the plots are 
given as maximum, most likely, and minimum expected. 

In Section [V] we justify some of the approximations 
and assumptions of the previous sections. We compare 
our results with others from the literature. Our notions 
of continuous background and unresolvable background 
are compared with the ones of previous work. We also 
show that the overlap function, which turns out to be a 
generalization of the duty cycle, is a proper quantifier of 
the resolvability and continuity of the background. 

The main conclusions and results are summarized in 
Section |VI[ Those readers who are short of time, or pri- 
marily interested in the main results, are suggested to go 
directly to this section. 



II. CHARACTERIZATION OF THE 
BACKGROUND 

A. Notation 

All magnitudes are measured in the frame of the cos- 
mological fluid, since massive objects that are not subject 
to external forces come quickly to rest with respect to this 
frame. 

We use the index "e" (for emission) for quantities mea- 
sured close to the system at the time of the emission of 
the radiation. For example, f e (emitted frequency) is 
the frequency of a wave, measured soon after its emis- 
sion, before the expansion of the universe stretches its 
wavelength. Quantities measured here and now (which 
are called observed quantities) have no index. The fre- 
quency of the wave of the previous example, measured 
today, is thus denoted by /. 

We use the indices "low" , "upp" , "min" and "max" to 
refer to lower, upper, minimum and maximum, respec- 
tively. 



B. Terminology 

We now introduce some terminology to avoid confusion 
or ambiguity throughout the paper. 

By system we mean a certain configuration of physi- 
cal objects that is a source of gravitational waves. An 
example of system is a pair of neutron stars inspiralling 
towards each other. 

We use the term ensemble for the collection of systems, 
all having similar properties and behaviour, formed from 
the Big Bang until now. An example is the population 
of coalescing neutron star binaries in the universe. 

By signal we refer to the total gravitational radiation 
emitted by a system. One system emits only one sig- 
nal, that can range over a large frequency interval and 
exist over a long interval of time. Despite the interval 
of time it lasts, a signal is assumed to be characterized 
by a certain redshift, which remains the same from the 
beginning until the end of the signal (in Section V B we 



discuss the validity of this assumption) . A signal is com- 
posed of signal elements, each characterized by a certain 
infinitesimal frequency interval. 

The total (gravitational wave) background is the col- 
lection of all gravitational waves in the present universe. 
It can be divided into different parts, according to dif- 
ferent criteria. For example, primordial and contempo- 
rary parts, resolvable and unresolvable parts or contin- 
uous and discontinuous parts. One can also divide the 
total background into many different parts, each con- 
formed by the contribution of a certain ensemble. By 
extension, we use the word background when referring 
to both the total background and to its different parts. 
Hence, we can talk about the background of neutron star 
binaries, which is the collection of signals of the ensem- 
ble of neutron star binaries. The part of this background 
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that fulfills the condition of unresolvability would be the 
unresolvable background of neutron star binaries. 



1. Primordial versus contemporary background 

The primordial background [5J 134) is composed of grav- 
itational radiation emitted in the early universe, at very 
large redshifts. It is analogous to the background formed 
by the cosmic electromagnetic radiation [35] - In the 
case of the latter, the radiation was released (when pho- 
tons decoupled from matter) roughly a hundred thou- 
sand years after the Big Bang. On the other hand, the 
primordial gravitational radiation was produced in a tiny 
fraction of the first second of the universe [35]. In this 
background might be hidden waves from inflation and 
cosmic strings [37| [38] . 

The other part of the total background is still being 
produced at present and thus we refer to it as the contem- 
porary background. It is made by many different systems 
that formed in the past (at redshifts less than ~ 20, which 
is the largest redshift assumed for the systems we study) 
and can also form today. Examples of such systems are 
coalescing binaries, rapidly-rotating compact objects or 
core-collapse supernovae (some references were given in 
Section [TJ) . 

In certain frequency ranges one can get a clear view 
of primordial signals, whereas in others the contempo- 
rary signals dominate. The detection of the primordial 
background would be the most direct way to observe pro- 
cesses of the very early universe. But valuable informa- 
tion would also be gained from the detection of the con- 
temporary background, for example about binary forma- 
tion and coalescence rates. Furthermore, predictions of 
the contemporary background can set constraints on the 
frequency ranges where the primordial one could be de- 
tected. The contemporary background is the main topic 
of this paper. 

In the literature, primordial and contemporary back- 
grounds are often called cosmological and astrophysical, 
respectively. We do not use this words to avoid ambigu- 
ity, since sometimes both terminologies are used together, 
for example when talking about cosmological populations 
of astrophysical sources [39 , which might be confusing 
for non-specialized readers. 



2. Unresolvable versus resolvable background 

It is useful to classify the components of the back- 
ground depending on their resolvability. We now briefly 
comment on this concept; precise definitions of what we 
mean by resolvable and unresolvable backgrounds can be 
found in Section [II Fl 

Signals spend different intervals of time at different 
ranges of frequencies. In the case of binaries, they evolve 
much more rapidly at higher than at lower frequencies. 
At lower frequencies they will thus overlap (i.e., they 



will be observed at the same time) more frequently than 
at higher ones. A frequency bin of width A/, which is 
the frequency resolution allowed by the detector and by 
the data analysis method, will often be filled by one or 
more signals at low frequencies. On the other hand, a 
frequency bin at high frequencies will be filled by one 
or more signals only sporadically, since signals are very 
short. 

An unresolvable part of the background exists as soon 
as a frequency bin is constantly occupied by an average 
of one or more signals. At frequencies where such a back- 
ground dominates, the waveforms of the signals cannot be 
distinguished from each other. When a waveform is not 
resolvable, one cannot obtain information from it, such 
as the characteristics of the system that emitted that ra- 
diation. Moreover, such waveforms cannot be subtracted 
out from the data. 

The rest of the background is the resolvable part. The 
waveforms of this part can be distinguished from each 
other. One can thus obtain information about the system 
by studying the waveform of the emitted radiation. 

For some authors, what we call the unresolvable back- 
ground is defined as the stochastic background, and the 
remaining gravitational radiation is called the total gravi- 
tational wave signal [16] . This is a reasonable definition, 
but conflicts with what is often called stochastic back- 
ground by many other authors (for example in [40] and 
other papers cited in Section [I]). A more precise definition 
for stochastic background can be found in [T]. 



3. Continuous versus discontinuous background 



We now briefly comment on the concept of continuity 



of the background. In Section II G we give precise defini- 
tions of what we mean by continuous and discontinuous 
backgrounds. 

A continuous background exists in a frequency interval 
[/low j /upp] (that can be, for example, the frequency win- 
dow of a detector) as soon as this interval is constantly 
occupied by one or more signals. If in that interval there 
are gaps between signals, or the signals occur sporadi- 
cally, the background is discontinuous. 

The condition of continuity tends to that of unresolv- 
ability when / up p — / low tends to A/. If a background 
is not continuous in an interval of frequencies, it is not 
continuous either in a smaller interval. Therefore, only a 
continuous background can be unresolvable. 

We point out that the continuity of the background is 
not as relevant as the resolvability. However, we include 
it in the paper for two reasons: 

First, the continuity has been often used in the lit- 
erature (for example in the already mentioned papers 
[131 Q3]) to define the different regimes of the background. 
Once we know how to account for the continuity, we will 
realize that it is not the right tool to be used. Instead, 
the resolvability is the fundamental property of the back- 
ground. 
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Second, the continuity can be used to determine how 
often the background is observed. Suppose we want to 
detect a signal of some kind, but there is a background 
covering the signal. If the background is discontinuous in 
a frequency band, sometimes that signal can be clearly 
seen, without any background. On the other hand, if the 
background is continuous, we need to subtract it from 
the data in order to see that other signal. As we show 
in Figure [6j ET has no continuous background from bi- 
naries in its frequency window; BBO, on the contrary, 
has a continuous background of binaries crossing its fre- 
quency window, so the subtraction of background signals 
is necessary in order to detect other kinds of signals (this 
problem has been treated in [4]). 



C. Heuristic proof of the lack of confusion noise for 
ground-based detectors 

We now justify, in a simple heuristic way, that there is 
no continuous background (and therefore no unresolvable 
background either) from binary systems at frequencies 
larger than 10 Hz. 

A neutron star binary takes ~ 10 3 s to evolve from 
10 Hz to the coalescence (which can be proved by using 
Equation (62)). The most realistic coalescence rate for 
these binaries (see Table [[} is of - 10 5 yr 1 , in the whole 
observable universe. This implies — 0.003 coalescences 
per second. One could naively say that, on average, one 



would see — 10 J x 0.003 = 3 signals. But that would 
only be true if all binaries were close to us, at redshift 
— 0. The farthest binaries (close to redshift — 5) that we 
observe today at frequency — 10 Hz, emitted at — 10 x 
(1+z) = 60 Hz (using Equation (|3|) and needed just — 8s 
to coalesce. An interval of time of ~ 8 s at redshift — 5 is 
now observed as an interval of — 8 x (1 + z) = 48 s (using 
Equation (|5|). This implies that an average of — 48 x 
0.003 ps 0.14 signals are observed. The number of signals 
expected to be observed is thus a number between 3 and 
0.14, which, after doing the proper calculation, turns out 
to be smaller than 1. Hence, neutron star binaries do not 
produce a continuous background at frequencies higher 
than 10 Hz. 

Other binaries whose signals could produce a continu- 
ous background in the frequency window of ground-based 
detectors are those containing a stellar-mass black hole. 
But these binaries have a smaller coalescence rate and 
need less time to coalesce, from an initial frequency of 
10 Hz. The product coalescence ratex duration of signal 
would thus be even smaller. Therefore they do not pro- 
duce a continuous background either. 

At frequencies larger than 10 Hz, hence, there is no 
continuous background from binary systems. Between 1 
and 10 Hz one could have a continuous background, but 
it turns out to be well below the realistic sensitivity of a 
ground- based detector (see Figure [6]). 



D. Cosmological model 

1. Metric 

We assume a spatially flat, homogeneous and isotropic 
universe, described by a Friedmann-Robertson- Walker 
metric, 



ds 2 



z 2 dt 2 e +a 2 (t e ) [dr 2 +r 2 [d9 2 + sin 2 (0) d<j) 2 }] , (1) 



where c is the speed of light. The time coordinate t e 
is chosen to be, for convenience, the look-back time: it 
is at present and to ~13.7Gyr at the Big Bang. The 
usual look-forward time would be just t' = to — t e , with 
which the form of the metric would not change. The 
dimensionless cosmological scale factor, a(t e ), is chosen 
to be a(0) = 1 at present. The coordinates r, 9 and are 
called comoving coordinates, because they move with the 
cosmological fluid. For example, two objects at rest with 
respect to the fluid, at positions r\ and r 2 (and equal 
values of and <fi) , have a comoving distance r = r 2 — T\ . 
This comoving distance remains the same at every future 
time. However, the physical (proper) distance between 
them is r p hys(ie) = o(t e )r, and changes with time as 
the universe expands. Setting r — at the Earth, the 
coordinate r of a distant object is its comoving distance 
from us. 



2. Redshifting 



by 



The definition of the cosmological redshift z is given 



1 + z 



«(0) 
a(t e ) ' 



(2) 



where, as already said, a(0) = 1. This equation gives 
the value of the scale factor at the time of emission of 
a graviton (or a photon) that is today observed with a 
redshift z. 

We now see how the expansion of the universe affects 
frequencies and energies of gravitational waves, as well as 
infinitesimal intervals of time (a derivation can be found 
in Section 4.1.4 of [S]). A frequency f e of a wave emitted 
by a system at a redshift z corresponds to 



/ = 



fe 



1 



(3) 



at the present time. Since the energy of a graviton is 
proportional to its frequency, a graviton emitted with an 
energy E e is today observed with 



E 



1 + z 



(4) 



An infinitesimal lapse of time dt e (emitted interval of 
time) measured at redshift z is today observed as 



dt = dt e [l + z\. 



(5) 



5 



From ([3| , it follows that an infinitesimal frequency inter- 
val df e emitted at redshift z is today observed as 



df = 



dfe 

l + z' 



(6) 



Similarly, from Q, an infinitesimal energy interval dE e 
corresponds to 



dE 



dE e 
l + z 



(7) 



today. 



3. Volumes 



Some important quantities in this work are defined as 
densities, i.e. per unit volume (by which we mean the 
spatial volume). Because of the expansion, it is conve- 
nient to speak of two different volumes: physical and 
comoving volume. 

The element of physical volume dV at fixed time t e in 
the metric |l]) is given by a 3 (t e )r 2 siri(9) dr d6 dcf). We 
consider only sources uniformly distributed in the sky, so 
we can integrate for all angles 9 and cj), obtaining 



dV = 4ira 3 {t e )r 2 dr. 



(8) 



The factor a 3 (t e ) accounts for the expansion in the three 
space dimensions. 

The element of comoving volume dV c is defined by 
dV c = a~ 3 (t e )dV, which, using Q, gives 



dV c = 4:Trr 2 dr. 



(9) 



Suppose eight galaxies (that, at large scales, can be 
thought as point-like) are placed at the vertices of a 
cube. With the expansion, the galaxies separate from 
each other and the physical volume of the cube increases, 
but its comoving volume remains always the same. Since 
we are assuming that all massive objects are at rest with 
respect to the fluid, no system enters or leaves a certain 
comoving volume. For this reason it is straightforward 
to measure densities (for example, the number density of 
systems) per unit comoving volume. 

It is useful to write the element of comoving volume 
in terms of redshifts, instead of distances. For that, we 
have to find a way to transform infinitesimal intervals of 
comoving distance dr into infinitesimal intervals of red- 
shift dz. Suppose we observe today two gravitons, one 
emitted at redshift z and the other at z + dz. Since both 
reach us at the same time, and both travel at the same 
velocity c with respect to the cosmological fluid, the one 
with larger redshift was emitted at a time dt e before the 
other, and thus at a comoving distance dr further away 
from us than the other. The path of the gravitons, mov- 
ing in a radial direction (d9 = d<f> = 0), is obtained by 
setting ds 2 = in Q, which gives 



dr 



-dt P 



(10) 



To write dt e in terms of redshifts, we use the definition 
of the redshift. One can differentiate Equation ([2| with 
respect to i e , obtaining dz/dt e = —a{t e )/a 2 {t e ). Using 
([2]) again and the definition of the Hubble expansion rate, 
H(t e ) = — d(t e ) / a(t e ) (where the minus sign appears be- 
cause of the use of a look-back time), one obtains 



dt P 



[l + z]H{z) 



dz. 



(11) 



Here, the Hubble expansion rate has been written as a 
function of the redshift, instead of the time. The form 
of H (z) is derived further on in this section. Introducing 
( 11 1 in ( 10 1, we obtain a relationship between dr and dz, 



dr 



1 



a(t e ) [1 + z]H(z) 



dz = 



H{z) 



dz, 



(12) 



where the terms a(t e ) and [l + z] have canceled out, using 
([2]). Finally, inserting (12) in Q, the element of comov- 
ing volume becomes 



dV c = 4irr 2 (z) 



H(z) 



-dz. 



Here, r(z) is obtained by integrating (12), 



r(z) 



H{z) 



(13) 



(14) 



Gravitons emitted between redshift z and z + dz define 
a shell of comoving volume given by dV c ■ 

The Hubble expansion rate can be written as a function 
of the redshift. For that, we use the Friedmann equation 
(see Chapter 27 of gU), 



H 2 (t e ) 



IttG 



P(te) 



a 2 (t e ) 



A 
3' 



(15) 



where G and A are the gravitational and cosmological 
constants, respectively. This equation is obtained from 
the Einstein equation, imposing the metric (TJ) and the 
stress-energy tensor of a perfect fluid (see Chapter 5 of 
|42) ) . We assume a spatially flat universe, which means 
with zero curvature, k = 0. The term p(t e ) is obtained 
from the equation of a perfect fluid of density p and pres- 
sure p (which is also obtained from the Einstein equa- 
tion) , 



P - 3H(t e 



P(te) 



P(t e ) 



0. 



(16) 



We can solve this equation considering a universe dom- 
inated by non-relativistic matter (also called dust), 
p(t e ) = p m (t e ), with the equation of state p m — 0. One 



obtains p m (t £ 



oc a 



(t e ), which, using (|2j), becomes 



Pm (z) = p° m (l + z) 3 , (17) 
where p^ is the current value of the density of matter. 



One can also solve (16) using the equation of state of 
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relativistic matter (radiation), p r = pc 2 /3. But one can 
prove that the resulting density, p r {z) = p®(l + z) 4 , dom- 
inates in the Friedmann equation only at very large red- 
shifts. The redshift at which both densities, p m {z) and 
p r {z) equate is z eq ps 3 x 10 3 (from [13]). Considering the 
rcdshifts involved in this work (z < 20) we neglect p r {z) 
compared to p m (z). Inserting (17) in (|15[) and rewriting 



the latter in terms of the present value of the Hubble 
expansion rate, Ho = [74.2 ± 3.6] kms -1 Mpc~ x (from 

mi 



where 



Here, 



H(z) = H £(z), 



£{z) = ^n m [l + z]3 + Q A . 



ft. 



m 2 



and CIa 



A 

ml 



(18) 



(19) 



(20) 



are two dimensionless quantities called the density pa- 
rameters of matter and dark energy, respectively. The 
most recent values for the cosmological parameters ob- 
tained by the Wilkinson Microwave Anisotropy Probe 
after seven years of measurements are given in [43) . We 
adopt a density parameter of matter Q m = 0.27 and of 
dark energy Q\ — 0.73. For simplicity we do not consider 
any uncertainty in these values. 

For a better understanding of the relationship between 
volumes and redshifts, we can see Figure [l] where a Pen- 
rose diagram [45] for the metric ([T]) is shown. Each point 
of the diagram represents a two-sphere at a certain con- 
formal time. The (look-forward) conformal time is de- 
fined by 



ch] = — [1 + z]dt e . 
The coordinates of Figure [I] are defined by 



r' = arctan(?7 + r) — arctan(?? — r) 
rj = arctan(ry + r) + arctan(ry — r) 



(21) 



(22) 



Introducing (21) in ([I]), the path of a graviton fulfills 
dr = cdr\. In the diagram we use c = 1, so that r = n and 
thus r' = rj' for all null paths. All gravitons that reach 
us today have traveled along the null path shown (the 
straight solid line connecting z — and z max ). This path 
cuts the horizontal axis at the moment of the Big Bang, 
fixing the horizon of our observable universe. For each 
infinitesimal interval of time dr\ (that describes the differ- 
ence between the emission of two gravitons that reach us 
today) there is a corresponding interval dz, along the null 
path, which represents a shell of infinitesimal comoving 
volume dV c . 



E. Spectral function 

Under the assumptions (discussed in pQ) that the back- 
ground is Gaussian, stationary, isotropic and unpolar- 



observable universe at present 
lines of constant distance 
lines of constant time 




comoving distance dr 



FIG. 1. Penrose diagram of a universe described by the metric 
([l]). The straight black line crossing z = and z = z max 
contains all the gravitons that we observe today. 



ized, all the information about it is contained in a di- 
mensionless function called the spectral function, defined 
by 



Pln(/) £ln(/) 



Pc 



C 2 Pc 



(23) 



where c is the speed of light and / is the observed grav- 
itational wave frequency. The present critical density of 
the universe is 



mi 

8itG' 



(24) 



This is the density that closes a universe with zero cos- 
mological constant. This means that p c is the density 



that, inserted in Equation (15) (using 17a = 0), gives a 
zero curvature (fc = 0) at present (t e — 0). The func- 
tion £in(/) is defined in such a way that the total energy 
density of gravitational waves in the present universe is 



£!„(/) din/. 



In the literature, £i n (/) is often written as 



de 
din/' 



(25) 



(26) 



Some other authors [21] [46] prefer not to use this nota- 
tion, arguing that it may lead to misunderstandings (one 
could mistakenly think that the energy density is dif- 
ferentiated) . Regardless of the notation, Equation ( 25 ) 
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must be fulfilled, so that e\ n (f)d\n f is the energy per 
unit physical volume of gravitational waves between In / 
and ln/ + dln/. Thus, £l(f) is the fractional energy den- 
sity of gravitational radiation, per logarithmic frequency 
interval, in the present universe. 

We first describe a system as seen by an observer close 
to it at the time of emission. The energy released in 
gravitational radiation between logarithmic frequencies 
In f e and In f e + d In f e is 



dE e = P e (f e )dlnf e 



(27) 



This defines P e (f e ): the energy spectrum of a system at 
the time of emission. From (27) it follows that P e (/ e ) = 
dE e /d\nf e . 

The energy spectrum of a system at the time of emis- 
sion can be related to the energy spectrum today. The 
present energy dE radiated by that system, with a loga- 
rithmic observed frequency between In / and In / + dln /, 
is 



dE = P(f)d\nf, 



(28) 



which defines P(f) = dE/d\nf. Applying ^ and Q to 
Equation (p8|), one obtains 



dE e 
l + z 



P(f)dlnf e 



Comparing it to Equation (27), 

P(f) = (1 + zy^eife) = (1 + zr'Petftl + z] 



(29) 



(30) 



The function P e (f[\ + z]) is explicitly given for the case 
of a binary system in Equation ( 59 1 . 



We now calculate the energy spectrum per unit comov- 
ing volume of an ensemble. The number of systems per 
unit comoving volume during a time dt e is 



dn = h{z)dt e 



(31) 



Here, n(z) = dn/dt e is the signal comoving density rate 
(number of signals per unit emitted interval of time per 
unit comoving volume). The comoving energy density 
spectrum of an ensemble is 



P(.f) = / P(f)dn = / P(f)h(z)dt 



This means, the comoving energy density spectrum p(f) 
measured today is what in Equation (23) was called £i n : 



the present energy density of gravitational radiation per 
logarithmic frequency interval (of a certain ensemble). 
Using J23l, p3l, and p2|, 



flnCf) 
PcC 2 



I f°° dt 
— / P{f)n{z)-ydz. 
PcC 2 Jq dz 



(35) 



In this formula, only the term dt e /dz depends on the 
choice of cosmological model. 

A similar derivation of ( [35] ) (using different notation) 
can be found in [33] . In that paper, the formula for the 
spectral function, called f2 gw (/), is given in Equation (5). 
The terms N{z) and {1 + z]~ 1 {f r dE gw /df r ] \ /,.=/[!.+*] cor- 
responds, with our notation, to h(z)dt e / dz and P(f), 
respectively. 

In Equation ( |33[ ) one can clearly see the assumption of 
a homogeneous universe, which is implicitly imposed by 
the metric ([I]). At any position within a shell of width 
dz there is the same number of systems. In other words, 
h(z) is the same at every point on a line of constant time, 
in Figure [I] 

The spectral function of an ensemble can be expressed 
more conveniently. We write it in terms of the energy 
spectrum at the time of emission, P e (f[l + z\), for our 
particular cosmological model. Using (11) and (30), 



«(/) 



p c c 2 H 



P e (f[l + z])n(z) 
[l + z] 2 £(z) 



dz. 



(36) 



The spectral function of the total contemporary back- 
ground would be the sum of the spectral functions of all 
different types of ensembles. 

But Q(f) does not include all redshifts and frequen- 
cies, since h(z) and P e (f e ) have support only for z e 
] and f e e [/mm, /max], respectively. The max- 
imum frequency /max is the one above which no more 
gravitational waves are emitted. The minimum frequency 
/min is the one below which the contribution in gravi- 
tational waves is dismissed. For example, neutron star 
binaries started to form at a redshift z max ~ 5 (~12 Gyr 
ago), are still forming at present, so z min ~ 0, and emit in 
a range of frequencies from ^0.01 mHz to ~1 kHz (these 



(32) ranges are justified in Section III). These limits in red- 



shift and frequency must be taken into account in the 



Recall that P(f) in general depends on z, according to integral of (36) 



(30). The integrals in 



from the Big Bang unti 



32 ) contain all systems formed 
Thus the limits of the 



To understand how ( 36 ) changes with the introduction 



today. 

time integral are (today) and to (the beginning of the 
universe). We can now change variables to write the 
previous integral in terms of redshifts, using the paths of 



gravitons as explained in Section II D 



P(f) 



P(f)n(z)^dz. 



(33) 



Since we have chosen a(0) = 1, the comoving volume and 
the physical volume are identical at present. Therefore, 



of these limits, it is helpful to make a plot of redshifts 
versus frequencies. Each horizontal line of such a plot 
gives the range of possible frequencies of a signal at a 
certain redshift. If we plotted on the horizontal axis the 
emitted frequencies f e , the limits / m i n , /max, z m - ul , and 
z max would define a rectangle, containing all the points 
(f e ,z) where both n(z) and P e (f e ) have support. But 
representing redshifts versus observed frequencies /, one 
obtains the plot of Figure [2] (which is no longer a rectan- 
gle). The shaded area represents the support of 0(/). 
We insert two redshift functions Z\ ov/ (f) and z upp (f) 



(34) in the integration limits of ( 36 ) , in such a way that the 




1+zmax 1+zmax 1+z 



FIG. 2. Redshift versus observed frequency. The spectral 
function of the total background of a certain ensemble has 
support only within the shaded area. 



integral is non-zero only in the shaded area of Figure [2] 
This is achieved with 



2low(/) 



/<T 



/min 



1. 



l+z„ 



/min 



</<lfc 
</ 



(37) 



and 



z upp(f) 




f<l 



+z„ 



1 + Z„ 



Jil 



(38) 



l+Zn 



With these limits, only signals emitted with frequencies 
between / m j n and / max , at redshifts between z m - ln and 
z max , contribute to 



Changing the integration limits in (36), the spectral 
function becomes 



«(/) 



1 



p c c 2 H 



P e (f[l + z])n(z) 
[l + z] 2 £(z) 



dz. 



(39) 



This formula gives the spectral function of the total 
background produced by an ensemble, measured at the 
present time. In the next section we generalize this for- 
mula to account for the resolvability of the signals. 



In several papers, the integral in Equation (39) con- 



tains an extra [1 4- z}^ 1 factor (see [18], where the origin 
of this misleading factor is explained). Recall the defini- 
tion of P e (f[l + z}) in Equation (27), and notice that the 
integral in ( 39 1 is not equivalent to, for example, that in 
Equation (9jof [15] . where the extra [1 + z]^ 1 factor is 
introduced [37] . 

Notice that the energy spectrum P e {f e ) does not de- 
pend on time. We are thus adding the contribution of 
each system as if it were instantaneous (this means, as if 
it were a point in the Penrose diagram of Figure [j}. This 
is justified if the inspiral times are much smaller than the 
cosmic timescales, so the time a system needs to evolve 



to /max is much less than the Hub- 



VB 



from emitting at f n 
ble time, Hq 1 sal3 Gyr. In Section 
this assumption is not always fulfilled, but it turns out 
to be irrelevant in practice. 



we point out that 



F. Resolvability of the background 

In this section we introduce the overlap function, 
M(f,Af,z), that allows us to define and quantify the 
resolvability of the background. 

We first define some quantities that are necessary for 
the definitions of the different parts of the background. 
Let B(f,zi,z 2 ) be the collection of signal elements with 
observed frequencies between / and / + df and with red- 
shifts between z\ and z%. Let r e (/ e , A/ e ) be the interval 
of time (measured close to the system at the moment of 
emission) that a system at z spends emitting between f e 
and f e + A/ e . Written in terms of observed frequencies, 
this interval of time is T e (f,Af,z). We define N(z) in 
such a way that N(z)dz is the number of signals pro- 
duced per unit emitted interval of time between z and 
z + dz. Since n(z) is the number of signals per unit emit- 
ted interval of time per unit comoving volume at redshift 
z, N(z) is given by 



N(z) = h(z) 



dVc 
dz 



(40) 



The value of N(z) at a certain redshift z can be con- 
sidered an average over an interval of time that is much 
longer than a typical observation time, but much shorter 
than the Hubble time. For the sake of simplicity let 
us assume that we know precisely this function, and 
that it gives the exact number of signals produced per 
unit emitted interval of time. For instance, if we have 
J* N(z)dz = lhour^ 1 , one signal is assumed to be pro- 
duced between redshift and z exactly every hour. 

Let us illustrate the resolvability with the following ex- 
ample: One signal is produced every hour between z and 
z + dz, i.e. N(z)dz = lhour -1 . Each signal spends one 
hour between / and / + A/, i.e. r e (/, A/, z) — lhour. 
Thus, whenever we see that frequency bin, it will be oc- 
cupied by r e (f, A/, z) x N(z)dz — 1 signal produced be- 
tween z and z + dz. If, for the same N{z)dz, we con- 
sider a different range of frequencies, where r e (/, A/, z) — 
2 hours, we will always see in that frequency bin 2 over- 
lapping signals, which will not be distinguishable. We 
can perform a similar calculation, considering all red- 
shifts between z% and z 2 '- J* 2 i~ e (f, A/, z) x N{z)dz gives 
the number of signals between redshift z\ and z 2 that 
overlap in a frequency bin. If that number is larger than 
1, those signals cannot be resolved. This leads us to the 
definition of the overlap function. 

The overlap function is defined by 



M{f,Af,z) 



,Xf,Af,z')N(z')dz'. (41) 
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It thus gives the expected number of signals with red- 
shifts smaller than z and frequencies between / and 
/ + Af. For example, A"(/, Af, z) = 1 implies that, as 
soon as one signal leaves a frequency bin, another signal 
enters it, so the bin is constantly occupied by one signal. 
Hence, JV(f,Af,z) > 1 implies that signals overlap in 
a frequency bin. We can impose N{f, Af, z) = Ao and 
invert this equation with respect to the redshift z. The 
obtained function, 



Af- 1 =Af- 1 (f,Af,Af ) 



(42) 



is the redshift such that all signals between / and / + Af 
with redshifts smaller than Af~ 1 (f, A/, Ao) sum Ao- To 
obtain an overlap larger than Ao at a certain frequency /, 
one has to consider only sign als from redshifts larger than 
Af~ 1 (f, Af, Ao). In Section |III E 3| we give a formula for 



In Section III 
N~ 1 (f, Af, Ao) for an ensemble of binary systems. 
We now give some relevant definitions: 
The total background of an ensemble between fre- 
quency / and / + df is B(f,z low (f),z upp (f)). One can 
assign a spectral function to it, f2totai(/)- 

The total background can be divided into two parts: 
resolvable and unresolvable. If there exists a certain z x 
such that z\ ow (f) < z* < z upp (f) and A"(/, Af, z») = 1, 
the unresolvable part is B(f, z*, z upp (f)), and the resolv- 
able part is B(f,Zi ow (f),z*). If there is no such that 
Zlow(/) < z* < z upp (/) and N{f,Af,z*) = 1, the re- 
solvable part coincides with the total background and 
the unresolvable part is the empty set. One can as- 
sign a spectral function to the resolvable part of the 
background, n r esoivabie(/), and to the unresolvable part, 

^ ^unresolvable (/)• 

The resolvable part dominates at a frequency / when 
^resolvable (/) > ^unresolvable (/) ■ When this happens, 
even if there is an unresolvable background present, it is 
weak compared to the background of the closer (stronger) 
signals, and thus the latter can still be resolved. On the 
other hand, the unresolvable part dominates at a fre- 
quency / when n unreso lvable(/) > ^resolvable (/)• In this 

case, even if there are some close resolvable signals, they 
cannot be resolved in practice, since they are obscured 
by the superposition of many weak distant signals. 

In Section |VD other possible criteria for the resolv- 
ability of the background are commented on. 

In Figure [3] we give an illustrative example to under- 
stand the definitions of the different parts of the back- 
ground. There we plot the evolution in time of the ob- 
served frequency of many similar signals, like the ones 
produced by an ensemble of binaries. The horizontal 
axis range is an interval of time of the order of a typ- 
ical observation time. This axis is divided in small in- 
tervals At, which is the time resolution. The vertical 
axis can be considered the frequency window of a hy- 
pothetical detector, with such a low instrumental noise 
that allows us to observe signals emitted at very large 
redshifts. This axis is divided into small intervals Af, 
the frequency resolution. Darker pixels in the plot repre- 
sent stronger backgrounds, i.e. with larger spectral func- 



(A/) 



(A/) 
(A/) 




time 



FIG. 3. Observed frequency versus time. Each line represents 
the evolution of one signal (like the one produced by a binary) . 
Closer signals, as well as the superposition of many signals, 
are plotted darker than distant individual signals. Three fre- 
quency bins are distinguished: in (A/)i the total background 
is completely resolvable, in (A/) 2 there is an unresolvable 
part and a dominating resolvable part, and in (A/)3 there is 
a dominating unresolvable part and a resolvable part. 



tion. The bin (A/)i is in a range of frequencies where 
the total background is completely resolvable: all signals 
can be clearly distinguished from each other. In (Af)2, 
an unresolvable part starts to contribute, but close bi- 
naries can still be clearly distinguished from each other, 
since the resolvable part dominates. Finally, in (A/)3 
the unresolvable part of the background dominates over 
the resolvable one. One should keep in mind that this 
example does not accurately follow the definition of un- 
resolvability, since the spectral function does not account 
for individual signals. 

We now generalize the formula of the spectral function 
to account for the resolvability of the background. We 



solve the integral in (39) for the signals that fulfill the 



condition AT(f, Af, z) > Ao- For that, we can retain the 



same upper limit of the integral, 
lower one, replacing z\ OVJ {f) by 



,(/), and change the 



upp(/)? 

z(f,Af,AT ) = { AT-\f,Af,Af ), f Pt 

upp(/ )l 



./ ^ /p ,131111 
lin — f — /p.max • 
./pjinax ^ / 

We have introduced four limiting frequencies: / Pjmax 
(Jp.min) represents the maximum (minimum) frequency 
at which the unresolvable part of the background is 
present, and /d,max (/d,mm) represents the maximum 
(minimum) frequency at which the unresolvable part 
dominates over the resolvable. Using (43 1 we obtain the 



spectral function of an ensemble with more than Ao sig- 
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nals per frequency bin A/, 

1 /"ZuppC/) 

n(/,A/,jV )=— J— / Pe{,f[l + Z]) 



PcC H h{f,Af,U ) 

[l + z] 2 £"(z) 



(44) 



This is the main equation of the paper and a generaliza- 
tion of Equation ( 39 1 with which we can distinguish the 



different regimes of the background. 

The unresolvable background is fully characterized by 
the spectral function fi(/, A/, Af ). It is not easy to de- 
termine whether the assumptions mentioned at the be- 
ginning of Section [lI E| are always fulfilled for such a back- 
ground. But it is clear that valuable information is lost 
when using the spectral function to characterize a re- 
solvable background, where signals can be individually 
distinguished. 

The spectral function of the unresolvable part of the 
background is, according to the definitions given at the 
beginning of this section, 



^unresolvable (/) = ^(/, A/, 1), 



(45) 



where A/ can be chosen as the inverse of the observation 
time. On the other hand, the spectral function of the 
resolvable part is 



resolvable 



,(/) = fitotaiCf) - n 



unresolvable 



where 



n tota i(/) = o(/,A/,o). 



(47) 



Here, £l{f, A/, 0) coincides with the given in Equa- 
tion (39), and the value of A/ becomes irrelevant. 



Another picture that illustrates the distinct parts of 
the background is in Figure [4j This graph is the same as 
that in Figure [2j but also represents the redshift function 
z(/, A/, Ao) that defines the frontier between the resolv- 
able (light-shaded area) and unresolvable (dark-shaded) 
parts of the background. 

The mathematical definitions of the limiting frequen- 
cies can be understood by looking at the graph in Fig- 
ure HI The frequencies /d.min and /d,max are the ones 
at which the resolvable and the unresolvable parts have 
equal spectral function, so 

^ ^unresolvable (/a ,min/max) ^resolvable (/d,min/max) ; (^8) 

The frequencies / p , m in and / Pima x are the ones at which 
the function M~ 1 (f,Af,Af ) intersects z upp (/), so 

,rmn/max3 

Af,Af ) = z upp (/ p 

,min/maxj- 

(49) 



In Section [III E 4| we calculate these limiting frequencies 
for an ensemble of binary systems. 

So far we have distinguished the regimes of resolvabil- 
ity by using the frequency resolution A/, but not the 



time resolution At. In Section IV El we show how to re- 
define the overlap function to account for the time res- 
olution. In practice, the effect of introducing a realistic 
At in the calculations turns out to be irrelevant for our 
work. 

In Section |"V F 3| we show that the overlap function is 
a generalization of what in the literature is often called 
the duty cycle, D(z). The latter is proven to be a good 
quantifier of the unresolvability of the background only 
for very short signals (bursts). Furthermore, we use the 
name overlap function and not duty cycle, because the 
latter may be confusing: D(z) can be greater than unity, 
unlike the typical duty cycles used in electronics or in 
gravitational wave detectors. 



G. Continuity of the background 

The overlap function can be used to characterize not 
only the resolvability but also the continuity of the back- 
ground. We now give some definitions, similar to the 
ones given in the previous section: 

Given a frequency interval [/i ow , /upp] (that can be the 
frequency window of a detector), the total background 
<6(/, z\ ow (f), z npp (f)) of an ensemble between frequency 
/ and / + df (where / low < / < / upp ) can be divided 
into two parts: discontinuous and continuous. If there 
exists a certain z* such that z\ ow (f) < z* < z upp (/) 
,(/), (46) and J\f(f, min(/ max , / upp ) - /, z*) = 1, the continuous 
part is £>(/, z* , z upp (/)), and the discontinuous part is 
S(/> Z\ ovr (f) ) %*)■ If there is no z* such that z\ ow (f) < 



2* < Z 



UPP 



(/) and Af(f,Af,z*) = 1, the discontinu- 
ous part coincides with the total background and the 
continuous part is the empty set. One can assign a 
spectral function to the discontinuous part of the back- 
ground, £ldiscontinuous(/), and to the continuous part, 

^ ^continuous (/)• 

The definitions of resolvable and unresolvable back- 
grounds are valid both for signals which frequency in- 
creases in time, such as binaries, and for signals which 
frequency decreases (for which one could change Af by 
—A/ in the definitions). On the other hand, the given 
definitions of discontinuous/continuous backgrounds as- 
sume that the frequency increases in time. For signals 
with decreasing frequency, the condition of continuity 
would be Af(f, f - max(/ low , / min ), z*) > 1. 

In the following, when talking about the continuous 
background, we will assume /i ow = and / upp = oo. 
This implies that any part of the background that is not 
continuous in this circumstance is definitely discontinu- 
ous, for any other choice of /i ow and / upp - Besides, the 
unresolvable background is necessarily continuous. 

So, the spectral function of the continuous background 

is 

^continuous (/) — ^(/) /max /j 1)- (50) 

The spectral function of the discontinuous background is 

^discontinuous \f) — ^total(/) ^continuous (f) > (^-0 
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Resolvable part ; Unresolvable part 
dominates i dominates 



Resolvable part 
dominates 



••?(/• A/. A o) 



\ \ Unresolvable ipart 



ziow(/)\ 



£upp(/) 



•< v Resolvable pjart 



fp,r. 



fd,i 



/d,i 



FIG. 4. Redshift versus observed frequency. Each horizontal line contains the possible observed frequencies of a signal. The 
light-shaded (dark-shaded) area represents the resolvable (unresolvable) part of the background. The redshift functions 2iow(/), 
z(f, A/, Ao), and z upp (f) are shown with dashed, dotted, and solid lines, respectively. The frequencies / p , m in and / p , m ax delimit 
the interval where the unresolvable part is present. The frequencies /d, m in and /d.max delimit the interval where the unresolvable 
part dominates. 



where f2 to tai(/) * s ^ ne same of Equations (47 1 and (39). 



III. 



MODELS FOR THE ENSEMBLES 



Our work is focused on the contemporary background 
produced by coalescing binary systems. These are sys- 
tems composed of two objects that inspiral towards each 
other, producing gravitational waves with an increasing 
frequency until they coalesce. In order to emit gravita- 
tional waves significantly, they must be sufficiently mas- 
sive and/or compact. Each binary is assumed to be iso- 
lated and describing an orbit of zero eccentricity. Its 
components are assumed to be non-spinning. 

We sort the binary systems into two classes: stellar 
binaries and massive black hole binaries. 

By stellar binary we mean a system whose components 
have masses of the order of a solar mass (or tens of it). 
We consider those stellar binaries formed by two stellar- 
mass black holes (from now on we call this type of binary 



BH-BH), a stellar-mass black hole and a neutron star 
(BH-NS), two neutron stars (NS-NS), a neutron star and 
a white dwarf (NS-WD) or two white dwarfs (WD- WD). 
The majority of the star formation rates in the literature 
vanish at rcdshifts larger than 5 or 6 (see for example [4"81 - 
[54]). If we expect no star formation at higher redshifts, 
no coalescence from stellar binaries is expected either. 
We choose then a maximum redshift for binary coales- 
cences of z max = 5. The minimum redshift is chosen 
z m in — 0, since these binaries can also form at present. 

Massive black hole binaries (from now on, MBH-MBH) 
are systems believed to exist in the center of many galax- 
ies [551 156j . Their components have masses that range 
from ~ 1O 2 M to ~ 10 10 M Q . We consider four differ- 
ent models for massive black hole formation, presented in 
[57] : two of them (called SE/SC, for small seeds and effi- 
cient/chaotic accretion) with light-seed black holes pro- 
duced as remnants of Population III stars, and two with 
heavy-seed black holes formed from dynamical instabili- 
ties in the nuclei of protogalaxies (called LE/LC, for large 
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seeds and efficient / chaotic accretion). These formation 
models allow coalescences at redshifts reaching z w 20. 



A. Energy spectrum 

We assume that the energy lost by a system when emit- 
ting gravitational radiation between f e and f e + df e is of 
the form 



dE, 



e.sys 



= -K[fe} b df e , 



(52) 



for real constants k and b. This formula is valid for all 
systems considered in this work (binaries) and also other 
systems (see for example the emission model of magne- 
tars in |17j). In particular, for a binary system, 



1 



hc^M 5 } 1 / 3 and b: 



(53) 



Here we have introduced the chirp mass M. of the binary, 
defined by 



M = 



[mim 2 



13/5 



mi 



m 2 ] 1 / 5 



(54) 



where m\ and mi are the masses of the two components 
of the binary. 

This energy spectrum is obtained by assuming that 
the energy of the system, as well as the separation of 
the bodies, varies slowly with time. This is valid as long 
as the orbit is far wider than the last stable one (see 
Equation (67 1). In these circumstances, the system is 



well described by the Newtonian equations of motion of 
two point masses in a circular orbit. 



We thus derive Equation (52) for binary systems us- 
ing Newtonian mechanics. The energy of the system, in 
terms of the separation s between the stars of masses mi 
and 77i2 is 



E„ 



1 Gm\vni 

2 s ' 



(55) 



We reduce the two-body problem to one fictitious body, 
of mass equal to the reduced mass of the system, 
mim2/[mi + 7712], suffering the same force as each of the 
real bodies. Applying Newton's second law, 



Gm\mi 



7711TO2 
mi + m 2 1 



[27T/ ( 



i2 
e, orbit J 



(56) 



where /e.orbit is the orbital frequency, that is related to 
the frequency of the gravitational waves in the quadrupo- 
lar approximation by 



fe 2 f e , 



rbit- 



(57) 



Introducing it in Equation ( 56 ) one obtains a formula 



that relates the separation of the masses with the fre- 
quency of the gravitational waves, 



G[mi + m 2 ] 



1/3 



(58) 



Replacing (58) in (55) and differentiating, one finally ob- 



tains the energy spectrum ( |52| with the values of k and 
b given in ( 53 ) . A more detailed derivation of Equation 



( 52 1 can be found in Chapter 4 of [H] 



What we need is the energy spectrum of the gravi- 
tational radiation at the time of emission, in terms of 
observed frequencies, P e (/[1 + z ]) (recall the formula of 
n(f,Af,Af ) in Equation @). Using (J3j) , ([27]), and 



dE e 
d\nf e 

F 2/3 



dE, 



e.sys 



d In f e 



= .fcMfe]- 1/3 



k/^[i + z; 



2/3 



(59) 



The absolute value has been introduced because P e (f[l + 
z]) is a positive quantity; dE e is the amount of gravita- 
tional wave energy within a frequency interval, regardless 
of whether the energy of the system increases or decreases 
with the frequency. Equation ( 59 1 must be inserted in 

PI, 



fi(/,A/,A/- ) = / L s dz > 

p c c 2 H h {f , Af x ) [1 + z] 4 / J £(z) 

(60) 

to obtain the spectral function of an ensemble of binary 
systems. 



B. Interval of time per frequency bin 

To calculate the overlap function we need the interval 
of time that a system spends emitting in a frequency 
bin. The frequency of the radiation emitted by a binary 
evolves in (look-forward) time following the relation 



dtp 



<)(, 



.8/3 



GM 



5/3 



tll/3 



(61) 



the derivation of which can be found in Chapter 4 of 



Integrating (61 1 between f e and f e + Af e one obtains the 



interval of time that the signal spends in that frequency 
interval, 



Te(/ e ,A/e 

where 



/; 



■8/3 



[fe + A/ e ] 



-8/3 



5c 5 



256tt 8 / 3 [GM] 5 / 3 ' 



(62) 



(63) 



The interval of time r e (/ e , A/ e ) can be written in terms 
of observed frequencies, 

r e (f, A/, z) = 8 2 Q[f, A/] [1 + z}- 8 '\ (64) 

where we have defined 

Q(/,A/) = r 8 / 3 -[/ + A/]- 8 / 3 . (65) 



13 



The function r e (/, A/, z) gives the interval of time that 
a signal, produced at a redshift z, needs to change from 
an observed frequency / to / + A/. We remark that 
this is an absolute (positive) interval of time, and not a 
look-back time. 

In certain conditions Q(f, A/) can be simplified, by 
performing a Taylor expansion around A/ = 0, 



For NS-WD, we use the same criterion as for WD- 
WD, but considering that the radius of the neutron star 
is negligible with respect to the radius of the white dwarf. 
Hence, the maximum frequency is 



/•NS-WD 
' max 



'G[mi 



m 2 \ 



(71) 



Q(/,A/)«-A// 



-11/3 



(66) where r\ is the radius of the white dwarf, that can be 



for A/ < /. 



calculated with ( 70 ) 



C. Maximum frequencies 



The energy spectrum of binary systems (Equation 



(52)) is assumed to be zero outside a certain frequency 
range [/ mill ,/ mM ]. We now present our choices of / max 
for each type of system. We omit the index e to simplify 
the notation, but one should keep in mind that / max is 
an emitted frequency. 

For all binaries that do not contain a white dwarf, / max 
is reached when both stars are as close to each other as 
s\ so . This is the separation at the last stable orbit (see 
for example Box 25.6 in |42) ) . 



^lso — 



6G[mi 



m 2 \ 



(67) 



which is 3 times the Schwarzschild radius of each star. 



Using Equation ( 58 ), the frequency of the last stable orbit 
is 



lso 
max 



i 



6V6ir G[mi + m 2 ] ' 



(68) 



where we have used the index "lso" to distinguish this 
maximum frequency from the following ones. 

For WD- WD, since the radius of a white dwarf is much 
bigger than its Schwarzschild radius, we assume that the 
maximum frequency is reached when both stars touch 
each other. This happens when the separation between 
them is r\ + r 2 , the sum of their radii. This separation 



corresponds to a frequency (using again Equation (58l) 

1 / G[mi + ra 2 ] 



fi 



WD-WD 
max 



T2\ 



(69) 



For r\ and r 2 one can use 



0.0112i?, 



©1 





-2/3 




mi 




.WiCha. 







2/3 



(70) 



where mcha «1.44M is Chandrasckhar's mass. This 
formula gives the approximate radius of a white dwarf 
as a function of its mass rrn. It is obtained from Equa- 
tions (27) and (28) (where there is an extra factor M & 
on the right side) in [55] . 



D. Minimum frequencies 

There is a certain minimum frequency, / m in, such that 
the gravitational radiation emitted by a system below 
this frequency is disregarded, because other mechanisms 
of energy loss are more effective. It is difficult to find 
a precise description of these mechanisms for each type 
of system. We adopt a simple criterion, for all stellar 
binaries, that fixes the value of / m i n : the interval of 
time between the beginning of the inspiral phase (when 
the binary emits at frequency / m i n ) and the coalescence 
(when it emits at / max ) cannot be larger than a certain 
interval of time T max . This condition is equivalent to 
Te( /min, /max - /min) = T max . Using Equation ([62]), one 
obtains 



/min — 



/i 



-8/3 
max 



-i -3/8 



(72) 



For all systems considered in this work we can reasonably 
perform the approximation 



/min 



T, 



max 



-3/8 



(73) 



As we did with the maximum frequencies, we omit the 
index e to simplify the notation, but / m ; n is always an 
emitted frequency. The assumption of a maximum inspi- 
ral time is justified in Section [VB| The maximum inspiral 
times chosen are T max — 12 Gyr for stellar binaries and 
75 Myr for massive black hole binaries. These choices are 
now explained. 



1. Maximum inspiral time for stellar binaries 

For stellar binaries, T max = t 5 — 0, where i 5 is the look- 
back time at which the first stellar binaries coalesced (at 
z ps 5). Integrating (11), 



dt = 



1 



o [l + z]H(z) 



dz « 12 Gyr. 



(74) 



This choice of T max is somewhat arbitrary and even leads 
to an inconsistency: only binaries that coalesced at small 
redshifts could have that much time to evolve from an 
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initial frequency f m in until the coalescence. Moreover, 
at frequencies close to the minimum one, the approxi- 
mation of short inspiral times compared to the Hubble 
time, commented at the end of Section |IIE| is not valid 
anymore. In Section VB we justify our choice of T max , 
the exact value of which turns out to be unimportant in 
practice. 



2. Maximum inspiral time for massive black hole binaries 

The process that leads to two massive black holes co- 
alescing can be briefly summarized in three main phases 
|59| : dynamical friction, gravitational slingshot and grav- 
itational radiation. When two dark matter halos contain- 
ing black holes merge, the black holes suffer dynamical 
friction |60j with the environment and sink to the center, 
forming a wide binary (with large orbital period). At a 
certain distance the dynamical friction phase ceases to 
be effective. Then the binary can continue to shrink be- 
cause of three-body interactions with surrounding stars 
|61j . These stars are ejected from the center and subtract 
some energy from the binary in the process. This phase 
is called gravitational slingshot because of the ejection of 
stars. Eventually the dynamical friction plus the sling- 
shot phases shrink the orbit enough, so that the binary 
can continue inspiralling until a coalescence in a finite 
interval of time by only emitting gravitational radiation, 
which constitutes the third phase. Other possible evo- 
lutions involving interaction with surrounding gas have 
been investigated in the literature 62 64 . 

We impose that the minimum frequency is the one 
at which the gravitational slingshot phase ends and the 
gravitational radiation starts to dominate (see the discus- 
sion in Section VC). As we now show, this condition is 



reasonably well fulfilled by imposing the same maximum 
inspiral time T max = 75Myr for all masses. 

The frequency at which the slingshot and radiation 
phases overlap is the one at which the energy spectra of 
the two phases are equal, 



dEe 
dfe 



dEe 
dfe 



(75) 



The variation of the energy of the gravitational waves 
with their frequency, in any of the two phases, can be 
written as 



dE e 

dfe 



dE e dt e ds 

dt e ds dfe 



dE e 

dte 



ds 

dte 



-I -1 



ds 

dfe 



(76) 



Here, dt e is an interval of time and s is the separation of 
the two black holes, which is the semi-major axis of the 
ellipse described. Since the orbit is assumed circular, s 
corresponds to the radius of the orbit. The term dE e /dt e 
is the same in both phases. Also ds/df e has the same form 
in the two phases. Only the evolution of the semi-major 
axis in time, ds/dt e , is different. Thus, instead of finding 



the frequency that fulfills Equation ( 75 1 , we can obtain 
the separation s at which 



ds 

dte 



ds 

dte 



(77) 



and then calculate the corresponding frequency using 
@. 

Following [55] (or similarly [5U]), we write the evolu- 
tion in time of the semi-major axis of a binary in the 
two phases. In the gravitational radiation phase, this 
evolution fulfills 



ds 

dte 



64G 3 m 1 m 2 [mi + to 2 ] 



5c 5 s 3 



whereas in the gravitational slingshot phase, 



ds 

dte 



Ho-«s 2 
' 27rr? ' 



(78) 



(79) 



In the latter, T-L is the hardening rate, ct* is the velocity 
dispersion of the stars in the bulge of the galaxy, and r c 
is the core radius (see [511] for more details). We use the 
value of "H reached in the limit of a very hard binary, 
% sa 15. It is known that there is a correlation between 
o~* and the mass of the massive black hole ?tibh hidden 
in the bulge (see [55] and [53]). This relation (from the 
most recent fits, by [57]) is 



logi ° ("iff) = Cl + C2 logi ° 



200 km s" 



(80) 



with (ci,c 2 ) = (8.12±0.08, 4.24±0.41). From this equa- 
tion we obtain ct^tobh) and use ttibh = mi + m-2 to 
account for the two components of the binary. The core 
radius r c , in the limit of a very hard binary, grows during 
the gravitational slingshot phase until it reaches 



3G[mi + m 2 ] ^ ( Gm 2 \ 
4^| n {ac^J ' 



(81) 



where m 2 is the mass of the lighter black hole. 

We now calculate the separation s at which both 



phases overlap. Replacing (78) and (79) in (77) 



s 5 In' 



/ Gto 2 



V4^ 2 



727rG 5 TOim 2 [mi + m 2 ] 
5HcV 5 



(82) 



This equation can be numerically solved for each pair 
of equal masses mi = m 2 = m, obtaining the sepa- 
ration (let us call it sr) at which the gravitational ra- 
diation phase starts to dominate. Using ( |58[ ) one can 
calculate the frequency fn that corresponds to sr. It 
turns out that the obtained dependence of fn with m 
is very accurately fitted by fmin{m), defined in Equa- 
tion (73), using T max « 75Myr. This is a numerical 
coincidence that eases further calculations. The origin 
of this coincidence is the following: omitting the log- 
arithm on the left side of Equation (82), sr oc m/er* 
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while cr* oc to 1 / 4 - 24 . This leads to sr 



,0.764 



Ac- 
cording to Equation (|58j), / oc m 1/,2 s -3 / 2 and therefore 
fn oc m~ 646 . On the other hand, according to Equa- 



tion (73), /, 



5 / 8 = m 64 . Therefore, the depen- 
min with m are almost the same. As 



dences of fn and /, 
a conclusion, setting a maximum inspiral time of 75 Myr 
is (almost) equivalent to considering only waves emitted 
during the gravitational radiation phase. 



E. Calculations for stellar binaries 

1. Coalescence rate 

The signal comoving density rate h(z), that was de- 
fined in Equation (31), represents, in the case we study 



now, the number of binaries per unit emitted interval of 
time per unit comoving volume that coalesce at a redshift 
z. We can thus call it the coalescence rate or simply rate. 

To obtain n(z), one could choose a star formation rate 
from the literature (which is usually a function of the red- 
shift) and transform it into a coalescence rate, for which 
a coalescence probability distribution is necessary. This 
procedure is followed for example in [68] . In Section V A 



we show that the use of a constant coalescence rate is 
well justified, given the large uncertainties in the local 
coalescence rate. Therefore, to simplify calculations, we 
assume a rate of the form 



n{z) 



< z < z min 

^min — ^ — ^max 



(83) 



for a real constant R. The values of R for each ensemble 
are given in Table [I] 

Some of the coalescence rates in the literature are es- 
timated only within our galaxy. We need to extrapolate 
those coalescence rates, given per Milky Way equivalent 
galaxy, MWEG -1 , to the rest of the universe. One sim- 
ple way to translate galactic rates into rates per cubic 
megaparsec, Mpc -3 , is explained in Section 3 of (S3]. We 
use the same conversion factor of [70], which is referred 

to iTn, 



1 MWEG = 0.0116 Mpc" 



(84) 



A similar factor is given in Equation (4) of [72 . The 
conversion ( 84 ) assumes that the blue- light luminosity of 
the Milky Way is 1.7 x 10 10 Lb,®, where L B ,® is the blue 
luminosity of the Sun, while that of the close universe 
is 0.0198 x 10 10 Lb, © per cubic megaparsec. All these 
factors are very uncertain, as discussed, for example, in 
|73j . We assume no uncertainty in (84) but then round 



2. Spectral function 

We now rewrite f2(/, A/, Ap ) in a simple way. Intro- 
ducing the constant rate R in (60), 



n(/,A/,JVb) = <Si/ 2/3 [g(z upp (f))~g(z(f, A/,A/o))] 
Here, 

Rk 



6l ~ 2 IT 
p c c z H 

and g(z) is the solution of the integral 

g(z) = f [l + z]-^ 3 £- 1 (z)dz. 



(85) 
(86) 



(87) 



We solve this integral semi- analytically in Section|HI E 5| 



3. Overlap function 

We obtain an explicit formula for the overlap function 
of binary systems with a constant coalescence rate. In- 
troducing (13) and (64) in (41), 



Af(/,A/,z) = 



4?r 



<5 2 Qa,A/)[i + zr 8/3 

2 



R 



H Jo 



E-\z")dz" 



-£~\z>) 



dz' 



This intricate equation can be rewritten to obtain a sim- 
ple expression for the overlap function, 

Af(f, A/, z) = S 2 S 3 Q(f, Af)lg(z) - g(z low (f))}. (89) 

Here we have defined 



5, = 4vri? 



H 3 



and 

g(z)= f[l + z}-' ; 



|[i + z]- 8/3 J £~\z!)dz! 



(90) 



8~ 1 (z)dz. (91) 



This integral cannot be analytically solved. One can in- 
vert (89) with respect to the redshift, obtaining 



the coalescence rates to one significant figure. 



(92) 

In Section |IIIE5| we give a semi-analytical solution for 
Af(f,Af,z) and Af-^^Af,^). 
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BH-BH 


BH-NS 


NS-NS 


NS-WD 


WD- WD 


Minimum i?/[Myr _1 Mpc~ 3 ] 


1 x 1(T 4 


6 x 1CT 4 


i x i(r 2 


2 x 1(T 2 


2 x 10 1 


Most likely i?/[Myr -1 Mpc~ 3 ] 


5 x 10" 3 


3 x icr 2 


i 


4 x 10" 1 


1 x 10 2 


Maximum i?/[Myr _1 Mpc~ 3 ] 


3 x 10 _1 


1 


9 


9 


5 x 10 2 



TABLE I. Minimum, most likely and maximum coalescence rates assumed for each type of ensemble. The coalescence rates 
of BH-BH and BH-NS are taken from [73], where they refer to [75] and |76| . respectively. For NS-NS, the values are taken 
from [77] (our minimum and maximum values are the minimum and maximum ones allowed by the uncertainties). The rates 
of NS-WD and WD- WD are taken from Table 1 in |78| . In Section |IV| we consider also the recent coalescence rates of BH-BH 

-3 M ,„-l 



predicted in [75], of R = 0.361q'26 Mpc 3 Myr 1 . The values given in the literature per Milky Way equivalent galaxy are 



converted using (841. All coalescence rates are rounded to one significant figure. 



4- Limiting frequencies 



The limiting frequencies / p 



/d 



/d 



and 



/p,max are defined in Section II F For the systems we 
study, / Pimin and / d ,min turn out to be close to / min /[l + 
Zmax] , which is the minimum frequency at which the spec- 
tral function has support. For simplicity, we assume 



fn 



,/p.i 



On the other hand, the frequencies /d.max 



(93) 



and /, 



p,max 



must be calculated using Equations (48) and (49), re- 
spectively. 

We now show how to calculate / Pima x- Inserting Equa- 



tion (92) in (49) 



Q(/p,max,A/) = 



Ms 



<52<53[5( z upp(/p,max)) - g(ziow(/p,max))] 



(94) 

One can obtain / Pima x by solving this equation. However, 
one can use a more convenient formula for / Pimax , that we 
present now. All stellar binaries satisfy / m in/[l + z m i n ] < 
ip,max, so ziow(/ P ,max) = ^min- Adopting a frequency 



p,max 



IS 



resolution Af — lyr , the condition Af <C f 
fulfilled by all stellar binaries. We can thus use the ap- 
proximation of Equation (66) in (94), obtaining 



/p,n 



8AfS 2 S 3 [g(z )~g(z min )] 
3A^ 



3/11 



: -l -g(z mi „) 



~3A^ 



3/11 



s < i 

5 > 1 
(95) 



where the dimensionless parameter H is defined by 
- Q( l+'zllx ' A /)^^3 [g(z max ) -g(z min )} 



A/ 



p,max- 



(96) 



The 



If S < 1, we have a simple formula for /, 
condition S < 1 is fulfilled by all stellar binaries that do 
not contain a white dwarf. For NS-WD and WD- WD, 



5 > 1, and one has to solve Equation (95) numerically. 



Similarly, one can obtain a formula for the limiting 
frequency /d,max, using Equations (48), (85) and (66). 
We point out that /d, m ax is by definition smaller than 
/ P ,max- In addition, one can show that /d,max// P ,max 
cannot be smaller than a certain factor F, so 



' max < I 

/*p,max 



(97) 



The value of this factor is 



F = 



9 {9 1 (|[ffOzmax) +g{z miD )})) -g(z min ) 



3/11 



9(Zn 



g(Zmin) 



(98) 

For 2 max = 5 and z nlin = 0, one obtains F m 0.6. All 
stellar binaries except WD- WD fulfill that /d,max ~ 0.6 x 
/ P ,max- For WD- WD, /d, m ax and /p, m ax are almost equal; 
moreover, they are almost as large as / max - This means 
that the total background of WD- WD is almost entirely 
dominated by its unresolvable part. 

One should notice that fi(/, A/, 0) is equivalent to the 
old definition of the spectral function, f2(/), in Equation 
@. Setting A/ = 0, the function A/, 0) be- 

comes ziow(/) (using Equation (92)). Then, the limiting 
frequencies / p , min and f p<max become / m in/[l + and 
/max/[l + ^min], respectively (see Figures [3] and g). Using 
Equation (43), z(f, Af, 0) becomes identically zi ovr (f), 
and thus Equations ((39| and (44) become equivalent. 



5. Semi- analytical solutions 

In order to obtain a semi-analytical solution for 
Af, A/o), we need two functions, g(z) and g(z), that 
fit accurately the numerical solutions of the integrals in 
Equations |(87| and ([911. 



A possible choice of the functions g(z) and g(z) is 

g(z) = ai arctan Q4 (o 2 z" 3 ) (99) 

and 

g(z) = ai arctan" 4 (a 2 z" 3 ) , (100) 
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for certain parameters (01,02,03,04) and (01,02,03,04) 
that must be numerically calculated. The optimal pa- 
rameters between z min = and z max = 5 are 

(01,02,03,04) = (0.5604,1.235,1.0047,0.8364) (101) 

and 

(01,02,03,04) = (0.07024,0.8658,1.3236,1.511). (102) 
These values can be used for all ensembles of stellar bina- 



ries, since the integrals in Equations ( 87 1 and ( 91 ) depend 
only on cosmological parameters. 

The semi-analytical formula for the overlap function, 
using (100), becomes 



Af(f, A/, z) = S 2 S 3 Q{f, A/) oi arctan^ (o 2 / 3 ) 



We invert it with respect to the redshift, obtaining 



(103) 



Af- 1 (f,Af,N ) = 



tan 



a 2 



AA 

aiJ 2< J 3 Q(/,A/) 



i/o* \ 1 1/03 



(104) 



Introducing it in ( 43 1 we obtain a formula for 
z(f,Af,M ). 



Finally, using (99), the semi-analytical formula for the 



spectral function of binary systems is 

n(f, Af,Af ) =5 iai f/ 3 [arctan O4 (02^ p (/)) 

- arctan ai (a 2 ^(/, A/,7V ))] . (105) 



The redshift function z upp (f) is given in (38). The lim- 



iting frequency / Pjma x can be calculated as explained in 
Section HIE 4 using the semi-analytical formula (100) 
for g(z). 



6. Mass ranges 

We calculate the spectral function of an ensemble as- 
suming that all similar objects have equal mass. For ex- 
ample, in the ensemble of NS-WD, all neutron stars have 
equal mass ttzns and all white dwarfs have equal mass 
towd- For this reason, given a range of possible masses 
for an object, we should not consider values of masses 
too different from the mean one. 

For a neutron star, we assume a mass in the range 
1.3 < mNs/Afo — 1-7. This interval is taken from [50] . 
where the lower limit predicted is (0.878 - 1.284) Mq, 
and the upper limit, (1.699 - 2.663) M . We use the 
largest mass of the lower limit and the smallest mass of 
the upper limit and round all values to two significant 
figures. The most likely value is the average of both 
limits of the interval. So, our choice for the masses of 
neutron stars is 



NS j m NS i m NS 



mcd 



(1.3, 1.5, 1.7) Md 



(106) 



The mass distribution of white dwarfs of spectral type 
DA, according to [81], is described by a Gaussian dis- 
tribution with mean /i =0.606 M Q and standard devia- 
tion a =0.135 Mq. The distribution of white dwarfs of 
spectral type DB has \i =0.758 M and a =0.192 M . 
Since we do not make a distinction between DA and 
DB white dwarfs, we calculate the Gaussian distribution 
that best fits the average of both distributions, obtaining 
/i =0.663 Mq and a =0.177 M Q . Similar results can be 
obtained using, for example, the distributions given in 
[82] . We assume a minimum mass of fj, — a =0.49 Mq 
and a maximum one of fi + a =0.84 Mq. Thus, 

(m$g, m$g, m%5) = (0.49, 0.66, 0.84) Mq. (107) 

For stellar-mass black holes, we calculate the mean \i 
and standard deviation a of the list of masses given in 
Table 1 of [53], obtaining \x = 7.8 M and a = 3.7 Mq. 
We assume for the minimum mass /i — a =4.1 Mq and 
for the maximum one fj, + a =12 Mq. Again, the most 
likely value is the average of both. Similar results can be 
achieved with the masses of Table 1 of [84] . The masses 
we use are, therefore, 



(rn 



m mcd m max 
BH i m BH ' m BH 



(4.1, 7.8, 12) M ( 



0- 



(108) 



F. Calculations for massive black hole binaries 

The masses of MBH-MBH range several orders of mag- 
nitude. It is reasonable to expect a very different number 
of signals produced by binaries of chirp mass 10 2 Mq than 
by binaries of 10 10 Mq. To be consistent with the given 
definition of ensemble (a population of many systems 
with similar properties and behaviour) , MBH-MBH form 
a superensemble composed of many ensembles, each one 
characterized by an infinitesimal range of chirp masses. 

The coalescence rate now depends on the chirp mass 
and the redshift. Instead of n(z) one now has a signal 
comoving density rate of the form n(M,z). This gives 
the number of signals per unit emitted interval of time 
per unit comoving volume per unit chirp mass. We do 
not have an analytical formula for n(A4 , z) [5S] . 

The spectral function of the total background of the 
superensemble is 



fitotal(/) = f /3 



Mn PP (z',f) 



where 
h(M,z) 



h(M',z')dM'dz', 
(109) 



i[GirM] 
9c 2 H$ 



5/3 



h(M,z)[l+z\-^£- x {z). (110) 



One can notice that ( |109 ) is the same as (60), just chang- 
ing n{z) by fi(A4, z)dM and integrating over chirp mass. 
The functions A^i ow (2,/) and _M up p(z,/) give, at ev- 
ery frequency and redshift, the minimum and maximum 
chirp masses that can contribute, respectively. In other 
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words, the interval [Mi ow (z, f),M. npp (z, /)] contains the 
chirp masses of those binaries which, at redshift z, have 
minimum frequency f min < / and maximum frequency 
/max > /■ They are obtained by inverting / m i„ (Equa- 
tion (73)) and / ma x (Equation (68)), respectively, with 
respect to M. Hence, 



Mi ow (zJ) = 
and 

Mupp(zJ) 



5c 5 



2567r 8 / 3 G 5 / 3 T ir 



6 VG 2 6 / 5 7T G 



3/5 



[/[! + *]] 



-8/5 



(111) 



(112) 



In the last equation we have used that, if the two masses 
of the binary are equal, then m\ = m 2 = m — 2 1 / 5 .M. 

The overlap function of the total background of the 
superensemble is 



fz rM upp (z',f) 

Af(f,Af,z)= / Q(/,A/) 

J Zmin JMi ov {z'J) 

x I 2 (M',z')dM'dz', 



(113) 



where 



h(M,z) 



5c 5 









X 


47T 








.H J 



2567r 8 / 3 [G7W] 5 / 3 
n 2 

£- 1 {z')dz l 



-Jl + z] 



-8/3 



Ho 



£~\z) 



n(M, z) 
■ (H4) 



Equation (113) is the same as (88), just changing n{z) 
by n(A4, z)dA4 and integrating over chirp mass. 

Section 4 of Hi describes a discrepancy between 
a semi-analytical calculation of the unresolvable back- 
ground of MBH-MBH and a Monte Carlo simulation. 
The discrepancy occurs because the semi-analytical ap- 
proach does not take into account the discrete nature of 
the systems. To account for it, they change the range of 
masses considered in the semi-analytical calculation. We 
now proceed in a similar way, to calculate the unresolv- 
able part of the background. 

The average number of signals with frequency equal or 
larger than / and chirp mass equal or larger than Ai is 



Af(f,M)= / Q(/,/xnax(M')-/) 

JM Jzi aw (M',f) 

x h{M',z')dz'dM'. (115) 

We impose that a signal, emitted at frequency / by a bi- 
nary with chirp mass M, can contribute to the continu- 
ous or the unresolvable background only if Af(f, M) > 1. 
The largest chirp mass A4(f) that contributes at fre- 
quency / is obtained by solving N{f 1 M{f)) = 1. We 
calculate A4{f) numerically, use it as the upper limit 
of the integral over chirp mass in Equation (113), and 



equate Af{f, A/, z) to 1: 



M(f) 

Q(/,A/) 

M la „(z'J) 

x I 2 {M',z')dM'dz' = 1. 



(116) 



Inverting the result of this equation with respect to the 
redshift, one obtains A/" _1 (/, A/, 1). Signals with fre- 
quency / emitted by binaries with chirp masses in the 
range [A4i ow (z. /), M] form an unresolvable background 
if their redshifts are larger than A/" _1 (/, A/, 1). 

The spectral function of the unresolvable background 
is, therefore, 

2/ r max /""^ (/) 

^unresolvable =/ / / I\{M.,z) 

JJV- 1 (/.A/,1) JMioAz'J) 

x dM'dz'. (117) 



Similarly, one can solve 

•A/X/, /max - /, z) - 



M(f) 

QUJr 

M low (z',f) 

x I 2 {M\z')dM'dz' = 



1 



/) 

(118) 



and invert it with respect to the redshift, obtaining a 
function Af-^f, /max - /, 1). Replacing 7V _1 (/, A/, 1) 
by Af~ 1 {f, /max - /, 1) in Equation (fTl7), one gets 



9. 



continuous 



= / 2/3 



x dM'dz' , 



M(f) 



h{M',z') 
(119) 



which is the spectral function of the continuous back- 
ground. 



IV. RESULTS 

The main results of this paper are presented in Figures 
[5] [6j [7j and [8j In Figure [5] we show the spectral func- 
tion of the total background of each ensemble. In Figure 
|6j the spectral function is plotted only in those regions 
where the background is continuous. Finally, Figure [7J 
which is the most relevant plot of the three, shows the 
unresolvable background produced by the different en- 
sembles, assuming Mo = 1 and A/ = lyr -1 . In these 
three figures, the values of masses and coalescence rates 
are the most likely ones. 

One can clearly conclude, from Figure [7j that ground- 
based detectors operate (and will operate) in a fre- 
quency range clean of confusion noise from binary sys- 
tems. Without taking into account other possible sources 
of unresolvable background, this frequency range could 
be a good scenario for the detection of primordial back- 
grounds. 

In Figure [8] we have plotted the contribution of each 
ensemble separately. For each ensemble, there are three 
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FIG. 5. Spectral function of the total background versus observed frequency. The contributions of the different ensembles 
are calculated with the most likely values of masses and coalescence rates. No restrictions in the duration of the signals are 
assumed in this plot, which means that also very short and sporadic signals are taken into account. As discussed in the text, 
the spectral function in such circumstances should not be compared to the sensitivity curves of a detector. 



different curves of Q(f, A/, Ao): one maximum, one min- 
imum and one most likely, depending on the values of 
masses and rates. 

For stellar binaries, the most likely expectation of 
f2(/, A/, Ao) is obtained by using the most likely 
chirp mass and coalescence rate. The upper curve of 
f2(/, A/, Ao) is the upper envelope of all curves that are 
obtained using the maximum rate and sweeping over all 
possible values of chirp mass. Similarly, the lower curve 
is the lower envelope of all curves obtained with the min- 
imum rate and sweeping over all chirp masses. 

For massive black hole binaries, the most likely curve 
is the average of the spectral functions calculated, as ex- 



plained in Section III F for each of the four models con- 
sidered. The upper and lower curves are 10 and 1/10 
times the most likely, respectively. These uncertainties 
have not been precisely calculated. Given the lack of 
observational information about many of the parameters 
involved, any accurate calculation of the uncertainties 
would be arbitrary. More precise errors are calculated in 
[IB] , based on the results of different theoretical models. 



The ranges of uncertainty given in [16] are similar to the 
ones we propose. 

One sees in Figure [7] that the unresolvable background 
is clearly dominated by the contribution of WD- WD, be- 
low - lO" 1 Hz, and of MBH-MBH, below - 10~ 4 Hz. 

The contribution of galactic binaries, which is believed 
to produce confusion noise in the frequency window of 
LISA, has not been included in the plots, since it can- 
not be calculated using fi(/). The spectral function is 
calculated assuming signals that are distributed homo- 
geneously and isotropically in the universe. This means, 
is related to the average density of gravitational 
waves in the universe. But within the galaxy the density 
is larger than the average. Moreover, galactic binaries 
are distributed anisotropically along the galactic disc. If 
one uses the spectral function to plot the contribution 
of galactic binaries, one is claiming that the density of 
gravitational waves in the universe is as large as the one 
inside the galaxy. Some papers in the literature which 
deal with the confusion noise produced by galactic bi- 
naries are |87l 188] . The most important contribution to 




FIG. 6. Spectral function of the continuous background versus observed frequency. The contributions of the different ensembles 
are calculated with the most likely values of masses and coalescence rates. The sensitivity curves of LISAQ ET (from |86]h 
BBO (from [4]) and the complete Parkes PTA (from [16) 1 are plotted for comparison. 

a Obtained using http://www.srl.caltech.edu/~shane/sensitivity/MakeCurve.html with the standard parameters. 



this background is the one by WD- WD. According to 
[55] . galactic WD- WD produce a background about an 
order of magnitude larger than that of extragalactic ones. 

One has to be careful when interpreting Figure [5] That 
plot gives us information about the averaged total energy 
density of gravitational waves produced by each ensem- 
ble. The curve of the total background of NS-NS, for 
instance, enters the window of ET, but that does not 
mean that ET will see a constant noise curve like that. 
The signals of NS-NS are, in that frequency range, short 
signals, that will often (but not constantly) be detected 
with ET. The effective sensitivity curve of ET is thus 
not affected by NS-NS. On the other hand, an unresolv- 
able background of WD- WD with a rate larger than the 
most likely one would certainly affect the sensitivity of 
LISA. To avoid misunderstandings we have not plotted 
the sensitivity curves of any detectors together with the 
total background. In Figure[6]we have plotted sensitivity 
curves just to show that the background is discontinuous 
in the frequency band of ground-based detectors. 

At frequencies close to the last stable orbit, the Newto- 



nian spectrum that we have calculated may differ consid- 
erably from the real one, since the assumption of slow or- 



bits made in Section III A is no longer fulfilled. Thus, the 
exact shape of the spectral function at such frequencies is 
not accurate. However, the continuous and unresolvable 
parts of the background lie safely at lower frequencies. 



A new prediction on the coalescence rate of BH-BH 
was published |79| during the writing of this paper. The 
rate given in that paper, of R = 0.36±qH Mpc" 3 Myr" 1 , 
is much larger than the one in Table |TJ This high estimate 
is based on the observation of two binaries, both contain- 
ing a stellar-mass black hole and a Wolf-Rayet star. Such 
rates have also been predicted by simulations [90 , con- 
sidering low-metallicity galaxies. We show in Figure [9] 
the total, continuous and unresolvable backgrounds, re- 
spectively, that such a rate would produce, assuming the 
same mass ranges for black holes given in Section |III E 6| 
In Figure [9] we also show the upper and lower limits al- 
lowed by the new rate. 




V. DISCUSSION 



On the coalescence rate of stellar binaries 



In Section Till E 11 we have assumed a coalescence rate 
that has the same value during all cosmological epochs 
(see Equation (83)). We now justify that this assumption 



is reasonable, considering the uncertainties in the local 
coalescence rate. 

In [T4J the coalescence rates of BH-NS and NS-NS are 
calculated as a function of the redshift (see Figure 2 of 
that paper). The rates peak at around z « 1 and then 
decrease, becoming zero between redshift 5 and 6. The 
difference between the rate at the peak and the local rate 
(at z = 0) is a factor of « 2.1, for NS-NS, and « 1.7Jbr 



BH-NS. One can calculate g(z) (using our Equation (87)) 



tral function and the overlap function are proportional to 
g(z) and g(z), respectively, the overall difference between 
using a constant and a non-constant rate would also be 
less than a factor of 2. On the other hand, the value 
of the local rate has an uncertainty of several orders of 
magnitude (see Table |l| . We thus consider that a factor 
of 2 is negligible compared to a factor of (at least) 100. 
In addition, we can see in Figure [T6| that our estimate of 
the total background agrees with that of [13 , which was 
calculated using a non-constant coalescence rate. 

Assuming that the rates of other type of stellar binaries 
have a similar behaviour than those shown in |14) . we 
can conclude that the use of a constant rate is a good 
approximation for stellar binaries. 



and g(z) (Equation EFti) introducing in the integrals a B - ° n the minimum frequency of stellar binaries 



normalized non-constant rate like the ones of Figure 2 of 
[13] . The obtained functions g(z) and g(z) differ from the 
ones calculated with a constant rate by less than a factor 
of ~ 2. The value of this factor would not change signifi- 
cantly if one used other rate functions (as pointed out in 
[55] regarding the value of ((1 + z) -1 / 3 )). Since the spec- 



In Section |IIID| wc have defined / min as the gravita- 
tional wave frequency such that a binary, emitting at this 
frequency, needs an interval of time equal to T max , the 
maximum inspiral time, to reach coalescence. But one 
could in principle find, for each type of binary, a more 
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FIG. 8. Spectral function of the total, continuous and unresolvable backgrounds of the different ensembles, versus observed 
frequency. In each plot there are nine curves: three of them are calculated with the highest values of coalescence rates, three 
with the most likely, and three with the lowest values possible. Three of the curves represent the total background, three 
the continuous part, and three the unresolvable part. In the case of NS-WD, the total and continuous curves are almost 
indistinguishable. The same occurs for WD- WD with the total, continuous and unresolvable curves. 



precise criterion to define / m i n . 

One could define / m i n , for example, as a function of 
the velocity kick that the components of the binary ex- 
perience at formation. This velocity kick, which can be 
provoked by a non-symmetrical supernova explosion, can 



push one component of the binary with enough energy in 
a direction opposite to that of the other component and 
disrupt the binary. So / m i n could be the frequency at 
which the orbital velocity equals the velocity kick. With 
such a criterion, using realistic values of these kicks [{JlT - 
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FIG. 9. Spectral function of the total (top), continuous (middle) and unresolvable (bottom) background versus observed fre- 
quency. The contributions of the different ensembles are calculated with the most likely values of masses and coalescence rates of 
Tableli] except for the case of BH-BH. The high rate of BH-BH, taken from the recent paper [79] , is R = 0.36±oJg Mpc" 3 Myr" 1 . 
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FIG. 10. Spectral function of the total, continuous, and un- 
resolvable backgrounds of the ensemble of BH-BH, versus ob- 
served frequency. This plot is analogous to that in Figure [8] 
but using the rate from [JJ] of R = 0.3612 ie Mpc" 3 Myr"\ 
instead of the one in Table III 



193] . one obtains too long inspiral times, in some cases 
orders of magnitude longer than the age of the universe. 

Our choice of T max is in fact almost as long as the age of 
the universe. Therefore, only binaries that coalesced re- 
cently could have had that much time to evolve from their 
formation, as commented on in Section [ill D 1| However, 
when considering long inspiral periods, one takes into 
account part of the contribution from binaries that have 
not yet coalesced. 

In Section |II E| we point out that the formula of the 
spectral function assumes short inspiral times, so that 
each signal starts and finishes at approximately the same 
rcdshift. But each system needs « 12 Gyr to complete 
the process, and the expansion of the universe is indeed 
relevant during that interval of time. We now investigate 
the effect of this apparent inconsistency. 

Our rate R accounts for coalescences (and not for 
births) of binary systems. This means that we are count- 
ing systems that are emitting at frequencies close to / ma x, 
the frequency of the coalescence. What we may be count- 
ing wrong are systems emitting at low frequencies. 

Suppose a binary, very close to us, that started inspi- 
ralling « 12 Gyr ago and coalesces right now. We only 
see the high frequency part of the spectrum, which is not 
redshifted. The waves emitted at the beginning of the 
inspiral (at low frequencies, «12 Gyr ago) are now far 
from us. But an observer located that far away would 
observe those waves today highly redshifted. Our mis- 
take, assuming short inspiral times, is to claim that the 
distant observer measures that low frequency radiation 
without any redshift. So the spectral function should be 
more redshifted (and thus have lower amplitude) at low 
frequencies. 

We now estimate below which frequency this effects 
starts to be important. For that, we assume the follow- 
ing: we assign wrong redshifts as soon as the difference in 



redshift between birth and coalescence of a signal is larger 
than 1. In units of time (using Equation (fllj)), a differ- 
ence in redshift of 1 implies timescales larger than «7 Gyr 
at redshifts close to zero and larger than «0.4 Gyr at red- 
shifts close to 5. To be conservative, we assume that these 
effects arc important when inspiral times are larger than 
0.4 Gyr. The lifetime of a binary is larger than 0.4 Gyr if 
its minimum frequency is lower than sa 4 x 10~ 5 Hz for 
BH-BH, « 7 x 10~ 5 Hz for BH-NS, « 1 x 10" 4 Hz for NS- 
NS, and w 2 x 10~ 4 Hz for NS-WD or WD- WD. These 
frequencies are in a range where the spectra of all stellar 
binaries are covered under an unresolvable background 
of MBH-MBH. 

We thus conclude that the exact values of the minimum 
frequencies are not relevant in practice. Furthermore, 
the assumption of short inspiral times is not fulfilled for 
stellar binaries at frequencies close to the minimum, but 
this docs not affect the results. 



On the minimum frequency of massive black 
hole binaries 



The minimum frequency of each massive black hole 
binary, as explained in Section flH F[ is assumed to be the 
frequency fn at which the slingshot and radiation phases 
overlap. This means that we dismiss the gravitational 
waves radiated during the slingshot phase. 

It turns out that the introduction of the slingshot phase 
in the calculations has a very small effect (well within 
the uncertainty ranges) in the spectral function of the 
superensemble, at frequencies below ~ 10~ 9 Hz. The 
reason is the following: for each ensemble of masses be- 
tween M. and M. + dM., the effect of introducing the 
slingshot phase is noticeable only at frequencies below 
/min(A^) (the one calculated using Equation (73) with 
Tmax = 75 Myr) . But the main contribution of each en- 



semble to the superensemble is at high frequencies, where 
they have larger spectral functions (because of the / 2 / 3 
factor). In the superensemble, the only appreciable low- 
frequency contributions are those from ensembles with 
the largest masses and with non-zero coalescence rates. 
Therefore, the effect of introducing the slingshot phase in 
the superensemble is noticeable only at frequencies close 
to /min(-M), when M. is in the range of large masses (of 
~ 10 8 — 10 9 M Q ). These frequencies are smaller than 
~10- 9 Hz. 



D. On the condition of resolvability 



In Section |IIF| we state that signals between / and 
/ + df with redshifts larger than such that z\ ow {f) < 
(/) and Af(f, A/, z*) — 1 are unresolvable. We 



< z 



upp 



are hence imposing a one-bin-rule: we are not able to dis- 
tinguish signals if there are more than one per frequency 
bin. Other authors suggest other possible criteria, such 
as the three-bin-rule or the eight-bin-rule [53]. Accord- 
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FIG. 11. Spectral function of the resolvable and unresolvable 
parts of the background of NS-NS, versus observed frequency. 
The values of masses and coalescence rates adopted are the 
most likely ones. The resolvable and unresolvable parts calcu- 
lated with the one- and with the eight-bin rule are compared. 



ing to these criteria, the condition of unresolvability is 
reached when each three (or eight) frequency bins are 
occupied by at least one signal. We now comment on 
how using one of these criteria would change our results. 

Imposing an eight-bin-rule makes the condition of un- 
resolvability less restrictive: signals become unresolvable 
at higher frequencies than for the one-bin-rule. The re- 
sults would be almost unaffected in the case of WD- 
WD, since the curve of the unresolvable background is 
almost as large as that of the total background. The 
spectral function of the unresolvable background for the 
remaining stellar binaries would be slightly extended to 
higher frequencies. We can calculate the spectral func- 
tion with the eight-bin rule, just by changing A/ by 
8 x A/, so fi unrcsolvab io(/) = fi(/,8A/,A/u). In Figure 
[TT]we compare the spectral functions of the unresolvable 
background calculated with the one- and the eight-bin- 
rule, for the case of NS-NS with the most likely values of 
masses and coalescence rates. 

One can note that imposing an eight-bin-rule, instead 
of a one-bin-rule, has the same effect of assuming an ob- 
servation time of eight years, instead of one. The ex- 
pected observation time of LISA is indeed three years; for 
the PTA, longer observation times are feasible. So, for 
MBH-MBH, redoing the calculations with the eight-bin- 
rule is compensated with the use of longer observation 
times. As pointed out in [T5], the unresolvable back- 
ground changes by less than a factor of 2 for observation 
times between one and ten years. 

We now discuss another possible definition of unre- 



solvable background. In Section II F we say that, when 
the unresolvable part of the background dominates, there 
still exists a resolvable part. We could consider this re- 
solvable part as also unresolvable. For that, we could just 



z(f,Af,Af ) = 



/p,i 



Af-\f,Af,Af ), 
^low (/), id ,mm 

AT-H/.A/.jVo), /d.max 



f ^ /p,min 
in — f — /d,min 

</</ d .max 

ix — f — /p,max 
7p,max ^ f 

(120) 

On the left side of Figure [l2j a plot of redshifts versus 
observed frequencies (analogous to that in Figure [4| is 
shown, using the new definition of z(f, A/, Ao)- There 
we see that, between /d, m in and /d,max, there is no re- 
solvable background. On the right side of Figur e [12 we 
show the spectral function obtained by inserting (120) in 



(44 ), for the case of NS-NS with the most likely values of 
masses and coalescence rates. The difference between the 
spectral functions with the old and the new definitions 
of z(f, A/, Ao) is just a small peak at frequency /d,max- 
Since the results are almost unchanged, we prefer the 
definition of unresolvable background given in Section 
IIF because the definition of z(f, A/, A/p) given in Equa- 
tion ( 43 ) is simpler than that in ( |120[ ) . 



E. On the definition of the overlap function 



In Section [Il F| we mention that the time resolution, At, 
should be taken into account to define the resolvability, 
but in practice it is not important. We now explain how 
At should be introduced in the overlap function and why 
it is not necessary, in the circumstances considered in this 
work. 

Signals of binaries evolve more rapidly at higher fre- 
quencies, and spend therefore less time in a frequency 
bin. Eventually, at frequencies and redshifts larger than 
certain / and z, the interval of time r e (/, A/, z) can be- 
come smaller than the time resolution At. When this 
happens, all signals spend effectively an interval of time 
At in each frequency bin, and no less than that (since 
smaller intervals of time cannot be distinguished). Thus, 
if the signals are unresolvable at / and z, they stay unre- 
solvable for any larger values of frequency and redshift. 
Equation (41) could be generalized to take into account 
this effect, 



A r (/, A/, At, z) = f max(r e (/, A/, z'),At[l + z}' 1 ) 

(121) 



'W/) 

x h(z') < ^ r dz' . 
az' 



The factor [1 + is necessary to compare our time 
resolution (which is an observed interval of time) with 
the interval of time at emission r e (/, A/, z). 

The effect of introducing At affects our calculations 
only if T e (f,Af,z) < At[l + z}- 1 when Af > 1, i.e., if 
there is more than one coalescence every At[l + z\~ l . 
Taking the highest rate of Table [I] (which is the maximum 
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FIG. 12. Redshift versus observed frequency (left plot, analogous to that in Figure |4| and spectral function versus observed 
frequency (right plot, analogous to that in Figure [8] corresponding to NS-NS with the most likely values of masses and 
coalescence rates). The difference between these plots and the ones in Figures [4] and [8] is that here, as soon as the unresolvable 
part of the background dominates, all signals become unresolvable. 



rate of WD- WD), we see that 



Am 



dz' 



is greater than one for At greater than ~ 1/9 s. A rea- 
sonable choice for the time resolution is the inverse of 
the sampling rate of a detector, which, in the case of cur- 
rent ground-based detectors, is much smaller than 1/9 s. 
Therefore, the generalization ( |121[ ) is not necessary; the 
overlap function is well defined by (41). 



By artificially increasing the time resolution by several 
orders of magnitude, we see the effect that the overlap 
function of Equation (121) produces in z( f, Af , At,Afp) 
(which is obtained by inserting Equation (121) in (43)) 
and in Q(f,Af,At,Af a ) (inserting ( |l2lj ) irTQ). This 
effect is plotted in Figure [l3j There we see that, above a 
certain redshift and a certain frequency, all signals con- 
tribute to the unresolvable background. With this exam- 
ple we see that a large time resolution would lead to the 
existence of an unresolvable background in the frequency 
band of ground-based detectors. 



F. Comparison with previous work 

1. Unresolvable backgrounds 

In Figure [14] we show the unresolvable background pro- 
duced by the superensemble of MBH-MBH, and the sum 
of the unresolvable backgrounds of all stellar binaries 
(which is almost equal to the background made by only 
WD- WD). These curves are compared with other pre- 
dictions from the literature. The curve (a) is obtained 
from |16j . using its Equation (14) with the mean values 
of the parameters in (45), (46) and (47). That formula is 
given in terms of the characteristic amplitude, h c , which 



is related to the spectral function by 



(122) 



In terms of the strain amplitude Sh(f), the characteristic 
amplitude is 



hc(f) = y/lsjjj. 



(123) 



The curve (b) is taken from Figure 4 of [15], where h IIOB 
also represents a characteristic amplitude. Finally, the 
curve (c) is taken from Figure 16 of [IT]. In all cases we 
find a good agreement of our predictions with those from 
the mentioned papers. 

In Figure [15] the unresolvable background of MBH- 
MBH is shown, calculated with the four different models 
). The unresolvable backgrounds calcu- 
[T6] are also plotted for comparison. 



(see Section III 



lated in [J5] anc 



2. Background of neutron star binaries 

In Figure[l6]we see that our estimate of the total back- 
ground of NS-NS is in good agreement with the one in 
[F3] . The curve (a) in Figure [16] represents what in [13] 
is called shot noise (see Figure 2 of that paper). In that 
work, the existence of a continuous (and Gaussian) back- 
ground is also asserted; this corresponds to the curve (b) 
in Figure 16 If we compare (b) with either our continu- 
ous or our unresolvable curves, we find a big discrepancy. 

We now explain the origin of this discrepancy. In Sec- 
tion 



II C we pointed out that signals of equal observed 



frequency but different redshifts need different intervals 
of time to coalesce. However, in Section 3 of [13], all 
binaries are assumed to spend the same amount of time 
in a certain frequency interval, leading to the conclusion 
that the background is continuous at high frequencies. 
The same has been claimed in similar papers, for exam- 
ple 03125]. 
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FIG. 13. Redshift versus observed frequency (left plot, analogous to that in Figure |4j) and spectral function versus observed 
frequency (right plot, analogous to that in Figure [8] corresponding to NS-NS with the most likely values of masses and 
coalescence rates). These plots (unlike those in Figures [2] and [8) are calculated assuming an unrealistic time resolution of 
At — 600 s. With such a large time resolution, an unresolvable background would be present in the frequency band of ground- 
based detectors. 




FIG. 14. Spectral function of the unresolvable background of MBH-MBH (dotted line) and of all stellar binaries (dashed line), 
versus observed frequency. These curves are compared with previous predictions from the literature, which correspond to the 
unresolvable backgrounds of: (a) MBH-MBH from |16| . (b) MBH-MBH from |15| and (c) extragalactic stellar binaries from 

na. 
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FIG. 15. Spectral function of the unresolvable background 
of MBH-MBH, calculated with the four different models (LE, 
LC, SE and SC), versus observed frequency. For comparison 
we include the curves of the unresolvable backgrounds of [T5] 
and [16]. 
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FIG. 16. Spectral function of the total, continuous and un- 
resolvable backgrounds of the ensemble of NS-NS (with the 
most likely values of masses and coalescence rates) , versus ob- 
served frequency. We compare these curves with those given 
in [13]. Curves (a) and (b) correspond to what in that paper 
is called shot noise and Gaussian background, respectively. In 
the text we explain the origin of the discrepancy between (b) 
and our continuous or unresolvable backgrounds. 



In a later work, [T3] , the continuity of the background 
is calculated in a similar manner as we do; the redshift 
of the signals is properly taken into account to measure 
the interval of time that they spend in the frequency win- 
dow of the detector. But in this paper, the continuous 
background is treated as unresolvable, which is incorrect, 
as we justify now. Suppose there is a continuous back- 
ground of NS-NS in the frequency band of ET, such that 
there is an average of a few signals present in the band. 
Even if these few signals are observed at the same time, 
they do not overlap in the frequency domain; the sig- 
nals can still be distinguished in frequency, so they are 
resolvable. 



3. Overlap function versus duty cycle 

In the literature, the so-called duty cycle is often used. 
It is defined by 



D(z) 



T P h(z')—r-^dz' , 

K ' dz' 



(124) 



where r e is the duration of a signal in the detector win- 
dow. If one assumes that r e is constant, as for example 
in [T3J [nni IHB] , D{z) does not give any valuable informa- 
tion (this has just been commented on in the previous 
section). If r e = r e (/i, /2 — fi,z), i.e., if r e is the time 
that each signal of redshift z spends in the frequency 
window [/i,/2] of a certain detector, D(z) characterizes 
the continuity of the background. But the property of 
the background that is indeed relevant is the resolvabil- 
ity, which is measured by the overlap function, defined in 
Equation (fill. 



An overlap function like the one in Equation (41) is 



useful for quantifying the resolvability of long signals. 
Now suppose that there is an ensemble of systems that 
do not emit gravitational waves during a long period of 
time, but rather in a burst. For such systems one cannot 
obtain an accurate function r e (/, A/, z). In this case, the 
resolvability can be quantified using the overlap function, 
by changing r e (/, A/, z) to r e , the typical duration of a 
burst. This overlap function would then coincide with the 
duty cycle. In the case that r e could be smaller than the 
time resolution At, one should rather use the generalized 
overlap function of (1211. 



The overlap function is therefore a generalization of 
the duty cycle, that can be used for short or long signals. 



VI. SUMMARY AND CONCLUSIONS 

We have reviewed basic aspects of the gravitational 
wave background. We have derived a formula (Equation 
(39)) for the spectral function, £l(f), for an ensemble 
of many similar systems emitting gravitational radiation 
at different times and locations. This formula has been 
generalized to account for the duration of the signals and 
the observation time. With the generalized spectral func- 
tion, fi(/, A/, TVo) (in Equation (44)), one can distinguish 



between unresolvable and resolvable backgrounds (Equa- 
tions ( 45 1 and ( |46| , respectively) , and between continu- 
ous and discontinuous backgrounds (Equations ( 50 ) and 
(51 1, respectively). 

The resolvability is a fundamental property of the 
background. An unresolvable background (often called 
confusion noise or stochastic background) is fully charac- 
terized by Sl(f, A/, A/"o). A resolvable background is com- 
posed of signals whose waveforms can be distinguished 
and in some circumstances subtracted out of the data. 
Precise definitions of resolvable and unresolvable back- 
grounds can be found in Section [II F| Figure [4] illustrates 
the different regimes of the background. 
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The resolvability is characterized by the overlap func- 
tion, Af(f, A/, z), which gives the average number of sig- 
nals, with frequency / and redshifts smaller than z, per 
frequency bin A/ (the frequency resolution). A formula 
for the overlap function is given in Equation (41). In 
Section |"V F| we have shown that the overlap function is 
a generalization of the duty cycle. The latter has been 
often used in the literature to quantify the continuity and 
even the resolvability of the background, leading in some 
cases to incorrect results. 

The continuity is a secondary property of the back- 
ground, which just gives an idea of how often the signals 
are present in the frequency window of a detector. The 
overlap function can also be used to characterize the con- 
tinuity of the background, as explained in Section [II G| 

We have calculated the spectral functions of the back- 
grounds of stellar binaries (those containing white dwarfs, 
neutron stars or stellar-mass black holes) and of massive 
black hole binaries. In Table U we have summarized the 
values, taken from the literature, of the coalescence rates 
of each ensemble. The ranges of masses assumed for neu- 
tron stars, white dwarfs and stellar-mass black holes are 
in Equations ( 106 1 , ( 107 ) and (11081) , respectively. A semi- 



with their uncertainties. The total, continuous and unre- 
solvable backgrounds, using the rate of BH-BH recently 
predicted in [75], are plotted in Figures [9] and 10 



analytical solution of the generalized spectral function 
has been derived for stellar binaries (Equation (105)). 



The calculations involving massive black hole binaries 
have been performed numerically, using the coalescence 
rates obtained with the four models of [57] , 

The spectral functions of the backgrounds produced 
by the different ensembles are plotted in Section [TV] over 
the frequencies of all present and planned detectors. The 
total, continuous and unresolvable backgrounds are plot- 
ted in Figures [5j [6] and [JJ respectively, with the most 
likely values of masses and coalescence rates. In Figure [8] 
the same curves are plotted separately for each ensemble, 



The unresolvable background is dominated by white 
dwarf binaries, below ~ 10 _1 Hz, and by massive black 
hole binaries, below ~ 10 -4 Hz. These backgrounds 
could enter the frequency window of LISA, PPTA and 
BBO. The continuous background of BH-BH, using the 
recent coalescence rate predicted in [75], becomes more 
important than the one made by NS-NS, especially in the 
band of BBO. The confusion noise produced by galactic 
binaries has not been shown in the figures, since it cannot 
be calculated using the spectral function. Some papers 
in the literature which cover this issue are I87H891 . 



Finally, with Figures [6] and [7J we conclude that present 
and planned ground-based detectors are in a frequency 
range where no continuous or unresolvable backgrounds 
from binary systems are present. Therefore, without con- 
sidering other possible sources of confusion noise, this 
band could be suitable for searching for primordial back- 
grounds. 
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